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ABSTRACT 

We report about the detection of 10 clusters of galaxies in the ongoing 
Swift/BAT all-sky survey. This sample, which comprises mostly merging clusters, 
was serendipitously detected in the 15-55 keV band. We use the BAT sample to 
investigate the presence of excess hard X-rays above the thermal emission. The 
BAT clusters do not show significant (e.g. >2 a) non-thermal hard X-ray emis- 
sion. The only exception is represented by Perseus whose high-energy emission 
is likely due to NGC 1275. Using XMM-Newton, Swift/XRT, Chandra and BAT 
data, we are able to produce upper limits on the Inverse Compton (IC) emission 
mechanism which are in disagreement with most of the previously claimed hard 
X-ray excesses. The coupling of the X-ray upper limits on the IC mechanism to 
radio data shows that in some clusters the magnetic field might be larger than 
0.5 fiG. We also derive the first log N - log S and luminosity function distribution 
of galaxy clusters above 15 keV. 
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Introduction 



Galaxy clusters are potentially powerful observational probes of dark matter and dark 
energy. However, the use of clusters to measure cosmological parameters becomes accessi- 
ble only when astrophysical uncertainties are well understood and controlled. Indeed the 
non-thermal pressure due to Cosmic Rays (CRs), magnetic fields and turbulence, is a source 
of systematic bias w hen cluster masses are estimated using the assumption of hydrostati- 
cal equilibrium (e.g. lEnsslin et al.lll997l ). The detection of clusters' X-ray emission above 
~20 keV is a fundamental step towards the firm grasp of these processes. 

It is well understood that clusters of galaxies contain a large amount of hot gas, called 
intracluster medium (ICM), that comprises 10-15% of their total mass. Already the first 
X-ray observations indicated the presence of this optically thin plasma, characterized by an 
atomic density of about 10~ 4 -10~ 2 cm " 3 and temperatures of the order of 10 7 -10 8 K (e.g. 



Felt en et al. 


1966; 


Catura et al. 


1972) 



Sarazin 


1988; 


Petrosian 


2001) 



However, evidences gathered at different wavelengths point to the existence of a non- 
thermal c omponent. In particular the detections of an extended synchrotron radio emis- 
sion (e.g. 



Willson 



197(1 lHarris fc Milevl Il978l; iGiovannini et al.lll993l; iGiovannini fc Ferettil 



200CH ; iKempner fc Sarazin 



2001 



Thierbach et al.ll2003[) and, more recently, of a possibly non- 



thermal extreme ultravio let (EUV) excess (ILieu et al.lll996l;lBowver et alll999l ; lBonamente et al. 



200ll ; iDurret et al.ll2002l ) and soft excess (e.g 



Werner et al 



20071 ) suggest the existence of a 
non-thermal X-ray component originating from a population of relativistic electrons. This 
scenario is confirmed by the detecti on of non-thermal emission in the hard X-ray spec- 



tra of a few galaxy clusters (see e.g. iKaastra et al.l 120081 ; iRephaeli et al. 



2008, for a com- 



plete review) . Still its actual presence and origin remain controversia l (IRenaud et al.ll2006a 
Fusco-Femiano et al.l 120071 ; IWerner et al.ll2007l ; iLutovinov et al.l 120081 ) . 



A non-thermal component could arise from a population of point sources (e.g. AGN as in 



Katzlll976l ; iFabian et al.ll 19 761 ; iFujita et al.ll2007l ) or from inverse Compton (I C) scattering of 
cosmic microw ave background (CMB) photons by relativistic electrons (e.g. 



Sarazin 1999) 



1999 



Timokhin et al 



2004 



Rephaeli 


1979 


K (e.g. 


Sarazin 



Sarazin fe Kemp ncr 20 00]) and synchrotron emission from ultra-relativistic electrons 



Inoue et al.ll2005l ; lEckert et al.ll2007l ). If the origin of the high-energy 



emission is IC scattering, then the presence of a large population of relativistic electrons 
(Lorentz factor >> 1000) is required. This population could have been ac celerated in shocks 



of different origin. I ndeed it could be associated t o merger shocks (e .g. IFujita et al.l 12003 



Brunetti et al.l 120041 ) . dark matter bow shocks (e.g. iBykov et al.ll2000l ). ram-pressure strip 
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ping of infalling ga laxies (e.g. Ide Plaa et al.ll2006l ). jets, Active Galactic Nuclei outbursts 
( Fujita et alj 2007, iu the ca se of radio mini-halos such as in Perseus cluster), accretion 



shocks (e.g. 
low 1 Gyr). 
electrons 



Inoue et al.ll2005l ). Non-thermal electrons loose energy on short timescales (be- 



Therefore some models consider a continued supply of primary accelerated 
(i.e. via first order Fermi mechanism), while others assume a constant in-situ 
re-acceleration via CR collisions or second order Fermi mechanism. 



If clusters are a large reservoir of non-thermal particles, then they should emit at higher 
energies, up to the 7-rays. Indeed, if CRs acceleration is taking place at the sho ck fronts 
then 7-r ays can be produc e d via IC, non-thermal bremsstrahlung and 7r ° decay (e.g. lRephaeli 
19791 : bar fc Shavivl Il995l : iReimer et aD I2OO3L 120041 : iBlasi et alil2007h A statistical upper 
limit on the flux above > 100 MeV was obtained by IReimer et al. J2003h . analyzing the 
emission from 58 clusters observed with EGRET. 

The role of CRs in the for mation and evolution of clusters of galaxies has been much 
debated. IChurazov et al.l (120071 ) suggest that in massive galaxy clusters hydrostatic equilib- 
rium is satisfied reasonably well, as long as the source has not experienced a recent major 
merger. However, in non-relaxed clusters the non-thermal pressure due to CRs, magnetic 
fields and micro-turbulence can affect the mass estim ates based on hydrostatic equilibrium 
(e.g. iMiralda-Escude fc Babullll995l ; iNagai et al.ll2007l ). This would lead to a higher baryonic 



to total mass ratio. Knowing the importance of CRs, the mechanisms that heat the ICM and 
the frequenc y at which it is shocke d, is crucial for the upcoming X-ray and Sunyaev-Zeldovich 
surveys (see 



Ando k Nagail 120071 ). 



In this paper we report the Swift/BAT all-sky detection of 10 galaxy clusters in the 
15-55 keV band. This constitutes the first complete sample so far detected at these energies. 
We use this sample to investigate the role of non-thermal processes in clusters. The structure 
of the paper is the following. In § [2] we describe the Swift/BAT observations and discuss 



the properties of each individual cluster (§ 12. 2p . In § 13.11 we provide, for all the clusters, 
constraints on the non-thermal emission as well as an estimate of the clusters' magnetic fields 
(§ 13.31) . The cluster source count distribution and the luminosity function are derived in §HJ 
We discuss the results of our analysis in § [SJ while § M summarizes our findings. 

We adopt a Hubble constant of H = 70 h 70 km s -1 Mpc -1 , Qm = 0.3 and Qa = 0.7. 
Unless otherwise stated errors are quoted at the 90 % confidence level (CL) for one interesting 
parameter and sola r abu ndances are determined using the meteoritic values provided in 
Anders fc Grevessd (Il989l ). 
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2. The BAT X-ray Survey 



The Burst Aler t Telescope (BAT; iBarthelmy et al.ll2005l ). on board the Swift satellite 
( IGehrels et al.l 120041 ). represents a major improvement in sensitivity for imaging of the hard 
X-ray sky. BAT is a coded mask telescope with a wide field of view (FOV, 120° x 90° 
partially coded) aperture sensitive in the 15-200 keV domain. BAT's main purpose is to 
locate Gamma- Ray Bursts (GRBs). While chasing new GRBs, BAT surveys the hard X-ray 
sky with an unprecedented sensitivity. Thanks to its wide FOV and its pointing strat- 
egy, BAT monitors continuously up to 80% of the sky every day. Results of the BAT survey 
( Markwardt et ap2005 : Ajello et al. 2008a ) show that BAT reaches a sensitivity of ~1 mCrab 
in 1 Ms of exposure. Given its sensitivity and the large exposure already accumulated in the 
whole sky, BAT poses itself as an excellent instrument for looking for the (faint) emission of 
galaxy clusters above 15keV. 



For the analysis presented here, we used all the available data taken from January 2005 
to March 2007. Since most of the cluster emission is expected to be thermal and thus rather 
soft, the survey chosen energy interval is 15-55 keV. The lower limit is dictated by the en- 
ergy threshold of the detectors. The upper limit was chosen as to avoid the presence of 
strong background line s which could worsen the overall sensitivity. The data screening was 
performed according to I Ajello et al.l (I2008al ). The all-sky image is obtained as the weighted 
average of all the shorter observations. The average exposure time in our image is 3 Ms, 
being 1.3 Ms and 5 Ms the minimum and maximum exposure times respectively. The final 
image shows a Gaussian normal noise. Source candidates were identified as excesses above 
the 5<T level. All these objects are then fit with the BAT point spread function (using the 
standard BAT tool batcelldetect) to derive the best source position. 



As shown in I Ajello et al.l (j2008al ) cros s-correlating the BAT sources with the ROSAT All- 
Sky Survey Bright Source Catalogue (jVoges et al.lll999l ) provides an easy and solid way to 
identify a large fraction (~70%) of them. Most of the uncorrelated sources are not present 
in the ROSAT survey because of absorption (either along the line of sight or intrinsic to the 
source). The unidentified sources a re targeted by Swift / XRT which in less th an 10 ks can 



pinpoint the exact counterpart (e.g. iTueller et al.l l2005al Jbl: iKennea et al.ll2005l ). 



The details about the complete source list will be given in an upcoming publication. Here 
we report about the detection of galaxy clusters above 5 a in the 15-55 keV band. 
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2.1. Clusters Identification 



Identifying clusters of galaxies as counterpart of BAT objects is not a straightforward 
process. Indeed, coded mask telescopes are rather insensitive to diffuse sources which extend 
over angles much larger than the projection of the mask tile on the sky (i.e. a few tens of 
arcmin for BAT). Even thoug h procedures exist to quantify the extent of d iffuse sources in 



coded mask instruments (see Renaud et al. 2006a b; Lutovinov et al. 2008, for the case of 



the Coma Cluster), their application is limited only to high signal-to-noise (S/N) objects. 
Given the extent of the BAT PSF (22' ), Coma is the only object whose emission is clearly 
extended in our investigation. Thus, for all other objects, the morphology of the source 
cannot be used to understand whether the BAT source is associated to the cluster or only 
to its brightest AGN. We therefore performed a spectral analysis (see § 12 . 2j) of those BAT 
sources that are spatially associated with galaxy clusters. All sources presented here show 
a significant thermal component that we interpret as thermal bremsstrahlung from the ICM 
and thus are securely associated with the proposed clusters. Our sample contains 10 galaxy 
clusters. Table [T] reports the position, significance, total exposure time and other details of 
all the detected clusters. 



For each g a 

m 



2.2. Spectral Analysis 

axy cl uster we extracted a 15-195 keV spectrum with the method described 



Ajello et al.l (j2008d ). Here we recall the main steps: the details can be found in the 
aforementioned paper. For a given source, we extract a spectrum from each observation 
where the source is in the field of view. These spectra are corrected for residual background 
contamination and for vignetting; the per-pointing spectra are then (weighted) averaged to 
produce the final source spectrum. Thus, the final spectrum represents the average source 
emission over the time-span considered here (2.5 years). The accuracy of these spectra is 
discussed in § 12 .31 

For all the clusters, we extracted a 0.3-10 keV spectrum using archival observations 
of XMM- Newton, Chandra, and Swift/XRT. Considering that for BAT all clusters, except 
Coma, are point-like objects we extracted (unless otherwise stated) all cluster photons within 
10' from the position of the BAT centroid. In most cases, this selection allows us to include 
most of the emission of the cluster. For those cases where there is clearly emission outside 
of our selection region^, we accounted for the missed flux using the beta profiles available 



In some cases the extent of the selection region is limited by the size of the CCD. 
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in literature. The details are given in the case-by-case section (§ 12. 4p . The level of the 
background was evaluated in regions of the CCDs not contaminated by the cl uster emission 



or using blank-sky observations (e.g lLumb et al.ll2002t iRead fc Ponmanll2003l ). In all cases 



we considered the systematic uncertainty connected to the background subtraction, in the 
0.3-10 keV band, to be 2%. All spectra were rebinned in order to have a minimum of 50 
counts (>7cr) per bin. 

As a standard procedure, we started fitting all the spectra with the most simple and 
plausible spectral model. In all cases this was a single-temperature thermal model with 
absorption fixed at the Galactic value. Only if the value of the x 2 /dof was greater than 1, 
we tried to add a second thermal model or a power law. In this case we chose the model 
which produced the best improvement in the fit (evaluated using the F-test) and the best 
residuals. 

Various authors have reported detections of hard X-ray excesses for some of the clusters 
present in our sample. For those cases where we do not detect directly such component, we 
tested whether our data are consistent with the reported non-thermal hard X-ray emission. 
This was done adding a power law to the thermal model used. We fixed the power-law 
index to 2.0, which is a value generally accepted for the non-thermal hard X-ray compo- 



nent generated by IC of relativistic electrons off CMB photons (e.g. iReimer et al.l 12004 ; 
Nevalainen et al.ll2004l ). We then let the power-l aw no r maliz ation vary until the A% 2 in- 
crement was larger than 2.7(6.64). According to lAvnil (119761 ). this gives the 90% (99%) 
confidence level on the parameter of interest. This allows us to investigate the level of 
non-thermal flux which is consistent with our data. 



2.3. Accuracy of BAT spectra 

When dealing with spectral features which are at the limiting sensitivity of a given 
instrument, it is important to make sure that all systematic uncertainties have been carefully 
taken care of. In order to test the reliability of our spectral extraction method, we extracted 
>160 spectra at random positions in the sky at least 30' away from the potential ( or detected) 



X-ray sources reported in the INTEGRAL reference catalog (lEbisawa et al.ll2003l ). The mean 
(raw) exposure of our spectral sample is 4.6 Ms. In each energy channel, the average flux 
is consistent within 1 a with zero as expected for pure noise and for efficient background 
subtraction. Moreover, the S/N distributions (i.e. flux divided by its error in a given energy 
channel; examples are shown in Fig. [1]) are all consistent with Normal Gaussian distributions. 
Both findings show that our spectra are trustable in the whole energy range (15-200 keV) 
and that uncertainties are well estimated. 
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Moreover, we can use the randomly extracted spectra to measure the average spectral 
sensitivity of BAT in a given energy channel. This is done deriving for each energy channel 
the standard deviation of the flux distribution. As shown in Fig. [21 the 3 a sensitivity in 
each energy channel is very close, except above 100 keV, to lmCrab. 



2.4. Individual Cluster Analysis 

2.4-1- Perseus 



Swift J0319. 8+4130 is certainly associated with the Perseus cluster (Abell 426). The 
BAT detection (see Fig. [3]) is well centered on the cluster. Perseus is one o f the most 
studied galaxy clusters an d its detection in X-rays date s back to the seventie s (IFritz et al. 
197ll ; iForman et al.lll972l ). XMM- Newton observations (jChurazov et al.ll2003l ) showed that 
the central region is contaminated by the emission of the AGN hosted by the brightest 
gala xy in Perseus, NGC 12 75. A hard X-ray com ponent has been detected with HEAO 
1 by iPrimini et al.l (119811 ). iNevalainen et al.l (120041 ) used BeppoSAX and previous RXTE 
measurements to prove that this non -thermal component is y ariable and must therefore be 
connected to the central bright AGN. ISanders fc Fabian! (120071 ) reported, using Chandra, the 



presence of no n-thermal X-ray emiss i on in the core of P erseus in correspondence of the radio 



mini-halo (i.e. iGisler &: Mileyl 1 19791 ; iGitti et al.l 120021 ). This non-thermal emission, which 



displays a power-law behavior with phot on index of 2.0, seems to exceed the flux of NGC 
1275 by a factor ~3 JSanders et alil2004h . 



The BAT spectrum shows evidences of an hard X-ray excess. Indeed, it can be fit 
by a steep power-law (photon index of 3.5±0.1 and Xred = 2-3) while it rejects a simple 
bremsstrahlung fit (Xred = 3-6). The fit improves ixled = 1-50) if we use a composite model, 
sum of the (bremsstrahlung like) gas emission and the (power-law like) AGN emission. The 
improvement of the fit is statistically significant as confirmed by the f-test probability of 
1.2xl0~ 2 . The best fit temperature is G^lfjjkeV and the photon index is 2. 5lj'n- If we 
fix the photon index at the value (1.65) determined by IChurazov et al.l (120031 ). we derive 
an extrapolated 0.5-8.0 keV luminosity of ~ 0.4 x 10 42 erg s _1 which is in agreement with 
the luminosity measured by XMM-Newton. This supports the idea that the hard-tail seen 
in the BAT spectrum is due to NGC 1275 and not to a non-thermal component originated 
in the ICM. Moreover, if we extrapolate, using a power-law with phot on index of 2.0, the 
non-thermal flux found in the 2-10 keV ranee by ISanders et al.l (120051 ) to the 50-100 keV 
band we get a value of 2.7x 10~ n erg cm -2 s _1 . This flux is a factor ~4 larger than the total 
cluste r flux observed by BAT in the same energy band. Recently iMolendi fc Gastaldello 
(120081 ) analyzed a long XMM-Newton observation and did not find evidence for non-thermal 
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emission. According to them the discrepancy between the Chandra and XMM-Newton 
results is due to a problem in the effective area calibration of Chandra. 

An XRT observation of 5.4 ks was carried out in July 2007. Given the size of the 
XRT CCD, we extracted all source photons within 6' from the BAT centroid. The surface 
brightness profile of Perseus is best describe d by the sum of a power-law and of a beta 
model. Adopting this model, as suggested by lEttori et al.l ( 119981 ). yields that ~94% of the 
total cluster emission fall s within our select ion. The joint XRT-BAT spectrum can be fitted 
by a sum of two APEC (jSmith et al.ll200ll ) models and a power law. The low-temperature 



component, which accounts for the cool core of the cluster, has a temperature of 3.0lo>keV 
and an abundance of 0.43^Q 2 g solar. The warmer component displays a temperature of 
6. 40ln 71 keV and an abun danc e of 0.3lln'|^ solar. T hese results are in line with the analyses 



of IChurazov et al.l (120031 ) and ISanders et al.l (120051 ). Both the power-law photon index of 
1.7^07 an d the the luminosity in the 0.5 -8.0 keV band of ~8x 10 _42 erg s -1 are compatible 
with the values found, for NGC 1275, by lChurazov et al.l (120031 ) and the ones determined in 
the next section. The photon index is slightly harder than the av erage photon index (2.0) of 



BAT AGN, however it is not unusual for radio-loud objects (e.g. Ajello et al. 2008cl ) 



The Nucleus of Perseus 

In order to study more in details the nuclear emission, we analyzed a 125 ks long XMM- 
Newton observation (observation 0305780101). We extracted the spectrum of the nucleus in 
a radius of 25" and evaluated the local background in an annulus around the source region. 
We note that the results presented here are not sensitive to the radius of the extraction 
region if this is in the 10"-30" range. The 0.2-9.0 keV spectrum of the nucleus is well fitted 
(x 2 /dof=960. 1/731) by an absorbed power-law model with absorption consistent with the 
Galactic one and a photon index of 1.60±0.02. Moreover, we find evidence (at the 95% 
CL) of a Ka Iron line with equivalent width of 90.2±45.0eV. An absorbed APEC model 
with a temperature of 12.6±0.7keV provides a worst fit (x 2 /dof=1167. 1/732) to the data. 
In particular, the absorption would be required to be lower than the Galactic one at 99 % 
CL. This fact, in conjunction with the presence of the Iron line, supports the evidence 
that the non-thermal emission in the nucleus of the Perseus cluster is produced by the 
central AGN. The non-thermal luminosities in the 0.5-8.0 keV and 2.0-10.0 keV bands are 
7.6lo2 x 10 42 er S s_1 an d 6.5^Q 2 x 10 42 erg s _1 respectively. In order to check these results we 
extracted a similar spectrum of the nucleus using Swift/XBT data and selecting an extraction 
region of 10" . The XRT data are compatible with the XMM-Newton one. Indeed, fixing 
the absorption at the Galactic value we find that the XRT data are compatible with a 
power-law model with a photon index of 1.6±0.1 and that the 2.0-10.0 keV luminosity is 



- 9- 



8.2_{ x 10 erg s _1 . Thus, the nucleus displays a moderate variability between the XMM- 
Newton and Swift/XBT observation epochs. This supports, once more, the interpretation 
that the non-thermal emission is produced by the central AGN. 



2.4.2. Abell 3266 



Swift J0431. 3-6126 is associated with Abell 3266. FigureHJshows that the BAT source 
is well centered on the cluster emission as seen by ROSAT. A 3266 (also known as Sersic 



40-6) was first detected in X -rays by the UHURU satellite (|Giacconi et al. 



1972 ). Accord- 



ingly to many authors (e.g. ISauvageot et al.l 120051 ; iFinoguenov et al.l 12006k and references 
therein) Abell 3266 recently underwent a major merger, p robably with a sub cluster tha t was 
stripped during the enc ounter with A 3266 dense core. iDe Grandi fc Molendil (119991 ) and 
Nevalainen et al.l (120041 ) observed Abell 3266 with BeppoSAX. The first group modeled the 
BeppoSAX broad-band spectrum (2-50 keV) with a simple optically thin thermal emission 
model at the temperature of 8.1±0.2keV, while the second group found a marginal evidence 
(0.8 a) of non-thermal X-ray excess. 



The BAT spectrum, shown in Fig. HI is consistent with the findings of lDe Grandi &: Molendi 



(]l999|). A bremsstrahlung model with a plasma temperature of 6.9^ 2 gkeV provides indeed 
a good fit to the data (x 2 /dof = 7.2/10). XMM-Newton observed Abell 3266 for 8.6 ks in 
September, 2000. The cluster is not centered on the EPIC-PN CCD. Thus we could extract 
only photons within a circular region of ~8' radius centered on the BAT centroid. In order 
to estimate the flux missed by our selection, we ado pt, for the cluster sur face brightness, 
a beta profile with (3=0.51 and core radius R c =3.1' (ISauvageot et al.l I2005T ) . According to 
our estimate 80 % of the total cluster flux is contained in our selection. Therefore, when 
fitting jointly the XMM-Newton and the BAT data, we use such cross-normalization factor. 
The combined XMM-Newton-BAT spectrum is well fitted by a single APEC model with a 
plasma temperature of S.O^o^keV and 0.41^^3 solar abundance. We derive a 99 % CL limit 
on the non-thermal 50-100 keV flux of 5.70xl0 _13 erg cm -2 s -1 . 



Extended radio emissio n correlated with A 3266 has been reported (IRobertson fc Roach 



1990; Brown & Burns 199 



adopt the radio data from 
flux density S 27QQ MHz 



). In order to es t imate the magnetic field (see § 13.11 Table 3), we 



Brown fc Burns! (Il99ll ). based on the Parkes catalogue, namely a 
1.070 Jy and a spectral index a = 0.95. 
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2.4-3. Abell 0754 



Swift J0908. 9-0938 is associated with the well studied cluster of galaxies Abell 0754. 



X-ray maps indic ate that A 075 
olent m erger fi.e. iHenrv &: Briel 



RXT E flValinia et al 



1999 



is far from hydrostatic equilibrium , experiencing a vi- 



1995 



Henriksen fc Markevitchl 119961). Its detection by 



Revnivtsev et al.l 12004 ) and BeppoSAX (jFusco-Femiano et al. 



20031 ) above 15keV make the association of the cluster with the BAT source secure. While 
the RXTE detections do not measure any significant hard X-ray excess, BeppoSAX detects 
a hard tail with a significant deviation from the thermal component in the 50-70 keV en- 
ergy range. It is worth noting that the BAT centroic|§ falls ~6' west of the brightest region 
of the cluster (see Fig. [5]). Chandra analysis of the gas temperature spatial distribu- 
tion shows indeed that the BAT position corresponds to regions of hot (TpsIO-15 keV) gas 
(jMarkevitch et al.l 120031 ) . The analysi s of XMM-Newton data confirms the existence of hot 
regions in the west part of the cluster (IHenry et al.ll2004l ). On the other hand, centroid shifts 



as a function of the waveband are a common indication of a merging cluster (jO'Hara et al. 
2004h . 



The BAT spectrum, shown in Fig. [5j is well fitted (x 2 /d°f — 6.3/9) by a single 



bremsstrahlung model with a plasma temperature of 9. 9^9 k keV. This is in g ood a gree 



ment with the temperatur e of 9A^^keV rep orted by iFusco-Femiano et al.l (120031 ) and 
9.0±0.13keV reported by IValinia et all Jl999h . The BeppoSAX 10-40 keV non-thermal 



flux of ~1.6xl0 12 erg cm 2 s 1 is consistent with the (90%) upper limit from BAT of 
6.5xl0 _12 erg cm" 2 s _1 . 

XMM-Newton observed A 0754 for 11 ks in May, 2001. The XMM-Newton-BAT data 
are well fitted by a single APEC model with a plasma temperature of S.Sto'jJkeV and 
0.29±0.03 solar abundance. Adding a power-law model, with photon index fixed to 2.0, 
improves the fit (F-test probability 4.6 xl0~ 9 ). The best fit temperature is 9.3±0.4keV and 
the non- thermal 50-100 keV flux is 7.6±^ x 10~ 13 erg cm s . The non-thermal flux in the 
10-40 keV band is 1.7l^xlQ- 12 erg cm s a nd is in goo d agreement with the non-th ermal 
flux measured by IFusco-Femiano et al.l (120031 ). However, IFusco-Femiano et al.l (120031 ) also 
discuss the possibility that the non-thermal flux be produced by the BL Lac object 26W20. 
This object lies ~24' away from the BAT centroid, and outside the XMM-Newton field of 
view, thus we can rule out that it is contributing to the detected non-thermal flux. 



However, we note that several point-like objects appear in the XMM-Newton image 



2 Fo r an 8 a detection, the expected maximum offset of the BAT centroid is ~2.5' (see Fig. 10 in lAiello et al 
2008ah . 
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and within 10' from the BAT centroid. A simple hardness ratio analysis reveals that the 
hardest object is located at R.A., Dec. = 09 09 13.7, -09 43 05.4. The likely counterpart 
is 2MASS 09091372-0943047 for which beside the magnitude (bmag = 20.0) nothing else 
is known. The XMM-Newton spectrum is extremely hard. It can be well represented, in 
the 0.1-10keV energy range, by an absorbed power-law with photon index of 1.23^ n ^ and 



-0.24 

x 10 21 atoms cm -2 . The source flux extrapolated to the 10-40 keV band 
" 1 . It is th us clear that this sing l e sou rce accounts for the 
non-thermal flux detected both by our and iFusco-Femiano et al.l (120031 ) analyses. 



absorption of 5.6jl|*| 
is (1.3±0.3)xl0~ 12 erg cm" 2 s~ 



Valinia et al.l (119991 ) and IFusco-Femiano et al.l (120031 ) derive a lower limit for the mag- 
netic field B of ~ 0.2 /iG and ~ 0.1 /xG respectively. Our es timate of B, reported in Tabled 
uses the VLA observations from lFusco-Femiano et al.l (120031) (<S , 13 65 MHz = ^6 mJy, a = 1.5 ) 



and is consistent with the results of iBacchi et al.l (120031 ) and of IFusco-Femiano et al.l (120031 ). 



2.4-4- Coma 



Swift J1259. 4+2757 is associated with the Coma cluster which is one of the best stud- 
ied cluster of galaxies. Coma (aka A 1656) is a particularly rich and symm etric merging clus- 



ter. It has been known as a diffuse X-ray and radio source since forty years (Felten et al 
Forman et aDl 1 972l : I Willsonll 1 9 70T) . The cluster hosts a powerful radi o halo dFeretti fc Giovanninil 



1966 



1998) and both BeppoSAX (IFusco-Femiano et al. Il999h and RXTE jRephaehl 
20021 ) revealed the existence of non-thermal hard X-ray emission. 



2001 



Rephaeli fc Gruber 



However, the detection of t 
positive BeppoS AX detections ( 



lis hard X-ray excess is sti l l quit e controversial. Indeed, the 



Fusco-Femiano et al.l Il999i . 120041) of hard X-ray excess wer e 



challenged by iRossetti fc Molendil (12004 l2007h . According to iRossetti fc Molendil (l2007h . 
the significance of the non-thermal excess changes (decreases) with the best-fit plasma tem- 
perature and only a certain set of assumptions (e.g. temperature o f the I CM) leads to a 
significant hard X-ray excess. However, recently IFusco-Femiano et al.l (120071 ). using different 
software analyses and studying a large set of background observations, were able to con- 
firm their previous finding. Independent ly of the BeppoSAX results, the RXTE detection 
(lRephaelill200ll ; iRephaeli fc Gruberl 120021 ) of the hard X-ray excess remains unchallenged. 



L ately, Coma has also been targeted by INTEGRAL (jEckert et al.ll2007l ; lLutovinov et al. 



20081 1. lEckert et all (120071 ) show, in their combined XMM-Newton-INTEGRAL analysis, the 
presence of a hotter region (gas temperature of 12±2keV as compared to 7.9±0.1keV of the 
center) in the south-west region. The authors favor the possibility that this emission is 
produced by IC scattering because its spatial distribution overlaps the halo of radio syn- 
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chrotron radiation. iLutovinov et all (120081 ) using INTEGRAL, ROSAT and RXTE data 
showed that the global Coma spectrum is well approximated by a thermal emission model 
only and found very marginal evidences (1.6 a) for hard X-ray excess. Thus, in light of these 
results the evidences for non-thermal emissions in Coma seem not conclusive. 

Coma is the only cluster in our sample whose extent is larger than the BAT PSF. The 
analysis of point-like sources in the vicinity of the Coma cluster shows that the PSF full 
width at half maximum (FWHM) is 22' while the FWHM of the Coma detection is 26' . 
Using a simple Gaussian profile for the surface brightness of Coma yiel ds a 1 a exten t in th e 
10'-15' range. This is in agreement with the morphological analysis of lEckert et al.l (120071 ). 
Moreover, from Figure the offset between the BAT and the ROSAT centroids is apparent. 
In deed, the BAT centr oid f alls ~4' west of the RO SAT surface brightness peak. As discussed 
bv lEckert et al.l (120071 ) and ILutovinov et all (120081 ) for INTEGRAL, the high-energy centroid 
coincides with a region of hot gas likely due to an infalling subcluster. 

Coded-mask detectors suppress the flux of diffuse sources and in order to recover the 
exact source flux and si gnificance one needs to develop dedicated methods for the analysis 
of extended objects (e.g. iRenaud et al.ll2006bl ). Given the fact that Co ma is the only cluste r 
'resolved' by BAT, a dedicated analysis will be left to a future paper (lAjello et al.ll2008dl ). 
However, we can extract the spectrum treating Coma as a point-like source. This translates 
into an analysis of the source emission within a radius of ~10' from the BAT centroid. The 
BAT spectrum is well fitted by a thermal model with gas temperature of 9.131J 3^ keV. 

XMM- Newton observed Coma several times. We analyzed an observation of 16 ks which 
took place in June, 2005. The XMM-Newton spectrum was extracted (as described in 
§ 12. 2p including all photons within 10' from the BAT centroid. Integrating the surface 
brightness profile deri ved by ROSAT (beta model with /3=0.74 and core radius R c =10.7'; 
Lutovinov et al.l 120081 ) shows that our selection includes ~75 % of the total Coma flux. A 
fit to the XMM-Newton-BAT spectrum with a single-temperature model does not yield sat- 
isfactory results (x 2 /dof = 1168.9/858). We then tried to add a power law to the APEC 
model. Adding a power-law model improves the fit (x 2 /dof = 905.5/856) and results into 
a well constrained photon index of 2.1llo'i3- However, this fit leaves evident ('snake'-like) 
residuals at low energy (see below for the residuals of all Coma fits). These residuals might 
highlight the presence of another thermal component. Indeed, we find that a satisfactory 
fit (x 2 /dof = 846.5/856) is achieved using two APEC models. The most intense compo- 
nent has a temperature of 8.4 01° 9^ keV and an abundance of 0.211° ^ consistent with what 
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found by lArnaud et al.l (1200 ll ) and ILutovinov et al.l (120081 ) . The low-temperature component 
(T=1.45^; 2 }keV and Z = O.O5(±O.O2)Z ) accounts very likely for one or more of the X-ray 
sources in the field of Coma. Indeed, an hardness ratio analysis of these X-ray sources shows 
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that t heir spectra are compati ble with thermal models with temperatur es in the 0.1-2keV 
range (jFinoguenov et al.l 12004 ). According to iFinoguenov et al.l (120041 ). these objects are 
(non-AGN) galaxies with a suppressed X-ray emission due to reduced star-formation activ- 
ity. Summarizing, we believe that the double-thermal model explains better the data with 
respect the thermal plus power-law model because: 1) it produces the largest improvement 
in the fit (i.e. largest Ax 2 ), 2) it better reproduces the low-energy part of the spectrum, and 
3) it accounts for all the point-like sources which are present in the XMM-Newton observa- 
tion. The best fit, sum of two APEC models, is shown in Fig. EJ The residuals to all the 
fits described in this section are reported in Fig. [7] while their parameters are summarized 
in Tab. [5j 

Our 99% CL upper limit in the 50-100 keV band is 1.70xl0 -12 erg cm -2 s -1 . However, 
we remark that this spectrum is representative only of the 10' radius region centered on 
the BAT centroid. Indeed, since the IC and the thermal emissions are propo rtional to the 



electro n density and to its square respectively (F/c oc n e and F thermal oc n 2 ; e.g. lSarazin et al. 



1998a]), their ratio (IC/thermal) is expected to increase with the distance from the cluster. 



Moreover, the lower density and larger sound speed (with respect the ph ysical conditions 



in the core) make CR acceleration more efficient in the cluster outskirts ( Pfrommer et al. 



20071 ). For these reasons and because Coma is an extended source for BAT, of which we 



analyze only the core, we cannot exclude the presence of a non-thermal component which 
arises in the outskirts of the cluster. 



2.4.5. Abell 3571 



Swift J1347. 7-3253 is likely as sociated with the Abell 3571 cluster, which has also been 
detected in the RXTE Slew-Survey ( jRevnivtsev et al.ll2004l ). Its symmetric morphology, see 
left panel of Fig. El an d temperature map indicate that A 3571 is a relaxed cluster (e.g. 
Markevitch et al.lll998l ). However, the radio structure, of the complex in whi ch A 3571 lays, 
suggests that this cluster is in the late stages of merging (IVenturi et al.ll2002[) . We note that 



Abell 3571 is known to have a moderate cool core ( Peres et al 



1998). Past and recent studies 



do not report evidences for non-thermal hard X-ray emission in Abell 3571. A fit to the BAT 
spectrum with a bremsstrahlung model yields a temperature of 6.9j1H ? keV (in agreement wi th 



the mean temperature of 6.7llo' 4 2keV measured with Chandra by lSanderson et al.l (120061 )) 



but the chi-square (x 2 ed =1.76) is relatively poor. The BAT spectrum shows positive residuals 
above 60keV which might reveal the presence of a hard tail (see Fig. [8]). However, given the 
low S/N of our spectrum, adding a power-law component does not impro ye the chi-square 
XMM-Newton observed Abell 3571 for 12 ks in July, 2007. According iNevalainen et al. 
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( 120011 ). the surface brightness of Abell 3571 follows a beta profile with /5=0.68 and core radius 
R c =3.85' . Therefore, our region of 10' radius includes approximately 93% of the cluster 
emission. This factor is taken into account when performing the joint fit of XMM-Newton 
and BAT data. The combined XMM-Newton-BAT spectrum, shown in the right panel of 
Fig. [8], is well fitted by an APEC model with a plasma temperature of 6.01±0.21 keV and an 
abundance of 0.34±0.06 solar. The total 2-10 keV flux of (8.0±0.3) x 10~ n erg cm" 2 s" 1 is 
i n good agreement with the value of 7.3±0.4 x 10~ n erg cm -2 s" 1 measured by BeppoSAX 



(INevalainen et al.ll200ll ). Even though statistically not required, a non-thermal power-law 
(photon index fixed to 2.0) is well constrained by our data. Indeed, we are able to derive a 
50-100 keV flux of (1.4±0.5) xl0~ 12 erg cm~ 2 s _1 . 

The radio flux density from the NRAO VLA Sky Survey is S 138Q MhLz = ^ 



(ICondon et al.lll998l ). We could not find any reference for the spectral index, so we adopted 
the value of a = 1.5 which leads to the lower limit listed in Table [31 We note that steeper 
spectrum gives a larger upper limit for the magnetic field (e.g. a = 2) would yield a lower 
limit twice as large as the previous. 



2.4.6. Abell 2029 



Swift J151 1.0+0544 is likely associated with the Abell 2029 cluster, which has also 



been detected at high energy by RXTE, Bepp oSAX and Chandra (jRevnivtsev et al.l 12004 



Molendi fc De Grandilll999l ; IClarke et al.ll2004L respectively). Left panel of Figure [9] shows 
that the BAT source is wel l centered on the cluster emission as seen by ROS A T. Abell 2029 



has a moderate cool core (ISarazin et al.lll998bl ; iMolendi fc De Grandilll999l ). IClarke et al. 



(120041 ) present an analysis of Chandra observations of the central region and find signs of 
interactions between the X-ray and the radio plasma. The unusual central radio source 
(PKS0745-191) morphology would be typical of a merging cluster. They suggest that A2029 
is a cluster that started very recently to cool to lower temperatures. 

The BAT data alone are well fit (x 2 /dof = 6.89/10) by a simple bremsstrahlung model 
with a temperature of lO.Gll'IkeV. An 8ks long XRT observation took place in September, 
2005. Given the extent of the XRT CCD, we extracted all the photons within a 6' from 
the BAT centroid. The surface brightne ss profile follows a beta model with /3=0.64 and 
core radius R c =1.8' (ISarazin et al.l Il998al ) . Integrating the beta profile up to 6' yields that 
95 % of the total cluster emission is included in our selection. However, for the case of 
Abell 2029 t he beta profile fails t o explain the inner 1.8' region which is characterized by a 
bright core (ISarazin et al.lll998al ). Thus, our selection might include an higher fraction of 
the total cluster emission. Indeed, BAT and XRT data are well fitted without the needs of a 
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cross-normalization constant. The BAT and XRT data are successfully fitted, by an APEC 
model with plasma temperature of 7.45+0.34 keV and 0.39±0.09 solar abundance which is 
consistent with the Chandra results (IClarke et al.ll2004j ). From the combined fit, we derive a 
99 % CL upper limit to the non-thermal flux in the 50-100 keV band of 1.27xl0 -12 erg cm -2 
s _1 . However, we note that the fit leaves positive residuals at high energy. We thus used a 
second APEC model, with abundance fixed at 0.4, to account for them. The F-test confirms 
that the second thermal component is detected at 99.85 % CL. The best fit temperatures 
are 9.61 2 ,' cJkeV and 4.1^'gkeV respectively. This fit is shown in Fig. [HI Abell 2029 has been 
targeted by ground- based TeV telescopes; however no TeV emission has been detected so far 
rtPerkins et al1l2006h . 



Condon et al.l (119981 ) found S 13S0 mjj z = 527.8 mJy. We adopted the value of a = 1.5 



which leads t o the lower limits on the magnetic field estimated in Table We note that 



Taylor et al.l (11994] ) obtained a lower limit on the magnetic field of 0.18/iG using observations 



of the central radio galaxy. 



2.4.7. Abell 2142 



Swift J1558. 5+2714 is associated with the Abell 2142 me rging; cluster. The de- 
tection in the 3-20 keV band by RXTE (jRevnivtsev et al.l 120041 ) makes the associ ation 



1998) and 



Markevitch et al 



of the BAT source w i th th e cluster rather strong. According to iPeres et al.l (J. 
De Grandi fc Molendi (ho02), Abe ll 2142 has a cool core that survived the merger. 
( 120001 ) and lSanderson et al.l (120061 ) . using Chandra observations, noted that the core of A2142 
has a complex structure, probably with a poor cluster enclosed in the halo of a hotter larger 
cluster. This would explain the lower temperature in the center, without the presence of 
a cool core. The left panel of Fig. [10] shows a point-like source located <4' from the clus- 
ter center. This object is the Seyfert 1 galaxy 2E 1556.4+2725. Given the distance, both 
objects, the cluster and the Syl, are not separated by BAT. 

The BAT data are well fit by a simple bremsstrahlung model (x 2 /dof = 7.96/10) with 
plasma temperature of lO.ljJ^keV. We analyzed an XMM- Newton observation of 800 s in 
conjunction with the BAT data. In this case, we extracted separately the spectrum of 
the cluster and the spectrum of the Syl 2E 1556.4+2725. The latter one shows an X-ray 
spectrum typical of a Syl object: i.e. absorption consistent with the Galactic one and 
photon index of 1.981;q; 14 . The extrapolated flux in the 15-55 keV range is 2.3xl0~ 12 erg 
cm -2 s 1 and it is well below the BAT sensitivity. Therefore we can consider negligible the 
Syl contribution in the BAT band. The surface brightness of Abell 2142 profile follows a 
beta model with /3=0.83 and core radius R c =4.2' ( IHenry & Briellll996l ). Integrating the beta 
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profile up to 10' yields that 97% of the total cluster emission is included in our selection. 
However, for the case of Abell 2142 the beta profil e underestimate the b rightness of the inner 
3' region which is characterized by a bright core (IHenry &: Briellll996l ). Thus, our selection 
might include an higher fraction of the total cluster emission. Indeed, BAT and XMM 
data are well fitted without the need of a cross-normalization constant. The cluster XMM- 
Newton-BAT spectrum is well fit by a simple APEC model with a plasma temperature of 
8.40lo 45 keV. The fit is shown in the right panel of Fig. [TUJ This is well in agreement with the 
temperatures of 8.8ln o keV and 9.0±0.3 keV measured by Chandra and GINGA respectively 



(IMarkevitch et al.ll2000l : Iwhite et al.lll994f ). From our fit the abundance is 0.27t H solar. 



Since no hard X-ray excess is detected, we report 99 % CL upper limits. Using a power-law 
with photon index of 2.0, we derive from the XMM-Newton-BAT data a 99% CL upper 
limit to the 50-100 keV non-thermal flux of 1.6xl0~ 12 erg cm -2 s _1 . The 99% CL limit on 
the non-thermal luminosity is 6.1xl0 43 erg s" 1 . T he marginal (~ 2 a) BeppoSAX detection 



of a non-thermal emission (jNevalainen et al.l 120041 ) is a factor 5 larger than our upper limit 
and thus incompatible with our data. 



The presence of a radio halo was already reported by lHarris et al.l (119771 ). iGiovannini fc Feretti 
(120001 ) measured S l4QQ = 18.3 mJy. In absence of a measured index a we adopt the 

arbitrary value of a = 1.5 to obtain the magnetic field constraint listed in Table [31 



2.4-8. Triangulum Australis 



Swift J1638. 8-6424, shown in the right panel of Fig. [TT1 is likely associated with the 
hot X-ray cluster of galaxies Triangulum Australis. This cluster at z =0.058 has a l ready 



been detected in the ROSAT, RXTE Slew and INTEGRAL surveys (IVoees et all Il999 
Revnivtsev et al.l 12004 ; IStephen et al.l 120061 ). In particular the detections by RXTE and 



INTEGRAL above 1 5 keV make this association certain . The Tria ngulum Australis cluster 
may host a cool core (jEdge et al.lll992l ; iPeres et al.lll998l ). However, iMarkevitch et al.l (119961 ) 
used the temperature and entropy maps from ASCA and ROSAT to find an indication of 
the probable presence of a subcluster merge r and argue that the cool gas in the core does 
not require a cooling flow. IMarkevitch et al.l (119981 ) found that a non-thermal component is 
more likely than a cooling flow. 

The BAT spectrum is well fitted (x 2 /dof = 5.68/9) by a simple brems strahlung model 
with p lasma temperature of 13. 4+37 keV. A similar temperature was found by lMarkevitch et al 
( 119961 ) in the centre of the cluster. 



XMM- Newton observed the Triangulum Australis cluster for 7480 s in February, 2001. 
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According to the beta profile reported by iMarkevitch et al.1 (119961 ). selecting photons within 
10' of the BAT centroid includes ~92% of the cluster emission. We thus employ such 
cross-normalization factor when fitting XMM-Newton and BAT data. The BAT and XMM- 
Newton data are consistent with a pure APEC model. From the best fit, shown in the 
right panel of Fig. [TTJ, we derive a plasma temperature of 9.30^30 keV and an abundance of 
0.301 

o!o7 solar. The XMM-Newton-BAT temp erature is in agreem ent wi th the mean values 



of 9.06+^ keV and 9.50±0.70keV reported bv llkebe et all (120021 ) and bv lChen et all (120070 



respectively. Using a power-law with photon index of 2.0, we derive a 99 % CL upper limit 
to non-thermal emission in the 50-100 keV band of Q.5xl0~ u 



erg cm 



2 s- 1 . 



Condon et al.1 (119931 ) report a 4.85 GHz radio source centered ~7' away from the BAT 



centroid. They find an upper limit of 33 Jy. We adopt this flux and the arbitrary value of 
a = 1.5 to obtain the magnetic field constraint listed in Table [31 



2.4-9. Ophiucus 



see Fig. [121) aw a y from one of the most studied 



Swift J1712.3-2319 lays only 1.7' ( 
galaxy clusters, Ophiucus, discovered by 
energies by BeppoSAX and INTEGRAL (INevalainen et al.1 l200ll; iBird et al.l I2006L r espec 



Johnston et al.l (119811). The detection at high- 



tively) makes the association with the BAT source certain. IWatanabe et al.l (j2001al ) used 
ASCA to measure the X-ray brightness distribution and temperature map. Considering 
the similarities with the Coma cluster, they conclude that Ophiucus is not relaxed and has 
likely experienced a recent merger. The BAT-derived plasma temperature of 9.5t^keV is 
i n good agreement with the values of 9. 6+n^keV and 9.0 1°^ keV measured by BeppoSAX 



( INevalainen et al.ll200ll ) and by Suzak u (IFuiita et al. 



2008 



An hard X-r a y exc ess was detected by INevalainen et al.l ( 200ll ) at a 2 a level. Very recently, 
Eckert et al.l (120081 ) using INTEGRAL confirmed this hard X-ray emission at an higher con- 
fidence level (4-6.4 a). The imaging capabilities of the instruments on-board INTEGRAL 
allowed the authors to conclude that the observed excess over the thermal emission is not 
originating from point sources (such as obscured AGNs) and is therefore non-thermal. This 
excess is marginally consistent with BAT data. Indeed, from our data we derive a 90 % 

1 while the 
i 



upper limit to the non-thermal component (20-60 keV) of 7.2x10 12 erg cm 2 s 



reported non-thermal flux observed by INTEGRAL is (10.1±2.5) xlO 12 erg cm 2 s 

We analyzed an archival Chandra observation of ~50ks. The observation, which took 
place in October 2002, was performed using the ACIS-S. Given its extent, the Ophiucus 
cluster is not entirely contained in a single chip. We thus extracted only those photons in a 
region of radius of 2.1' around the BAT centroid. The region extent is dictated by the size of 
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the chip. When performing a simultaneous fit with BAT data, we must therefore account for 
the flux which falls outside of the ACIS-S chip. Assuming that the surface density follows 
a beta profile and adop ting the values of 3 =0.64 and c ore ra dius of R c =3.2' as found by 
Watanabe et all (l2001bl ) and confirmed by lEckert et al.l (120081 ) , we derive that only ~52 % 



of the total cluster flux is included in our selection. If we let the cross-normalization of the 
BAT and the Chandra data vary, we derive that the Chandra data show a normalization 
(with respect the BAT ones) of 53^g % which is in good agreement with th e 52 % derived 



above. Thus, we fix the cross-normalization factor at 52 %. Moreover, as in Blanton et al. 



(120031 ) we accounted for the uncertainty in the background subtraction adding a systematic 
uncertainty of 2 %. The joint Chandra-BAT spectrum is well fitted by a single APEC model 
with a temperature of 9.93^o;t4 keV and abundance of 0.52±0.03. Using a power-law with 
photon index of 2.0, we derive a 99 % CL upper limit to the non-thermal emission in the 
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50-100 keV and 20-60 keV bands of 2.8x10-" erg cm 
respectively. The INTEGRAL detection is inconsistent 



-12 — 2 —1 

erg cm s 



s 1 and 4.5x10 



'2a) with our upper limit. 



The Oph iucus cluster is assoc iated in the radio domain to the extended radio so urce 
MSH 17-203 (1 Johnston et al.lll98ll ). The most recent radio data date back to 1977 ( jSlee 
19771 ) and report S mQ ]\/[g z = 6.4 Jy and a = 2, which we use to produce the lower limit 
on the magnetic fie ld reported in Tab. [31 The results do not change if we use older radio 
measurements (e.g. iMills et al.lll960l ; I Jones fc Finlayill974l ; ISlee & Hi gginsl 1 1 9 75l ) . 



24.10. Abell2319 



Swift J1920. 9+4357 is certa inly associated with the massive Abell 2319 cluster, that 
is undergoing a major merger (e.g lO'Hara et all 120041 ) . The BAT centroid (see left panel 
of Fig. [T3|) lies ~2' north-west of the peak of the ROSAT emission. Indeed, Chandra ob- 
servations reveal at the same position a region of hot, ~12keV, ga s while at the position 
of the ROSAT peak there is likely a cool core (lO'Hara et all 120041) . Abell 2319 has been 



detec ted above 10 keV by BeppoSAX and RXTE (IMolendi et al 



19991 : iGruber fc Rephaeli 



20021 . respectively). These two measurements are symptomatic of the uncertainty related to 
the hard X-ray detection claims from non- imaging instrume n ts and the inherent uncertainty 
from source contamination. Indeed, IMolendi fc De Grandi Jl999h repo rt that no hard-tail 
emission is present in BeppoSAX data, while IGruber fc Rephaeli! (120021 ) find that a power- 
law component can explain some residual features in the 15-30 keV energy range. The BAT 
data favors the thermal scenario. Indeed the best fit to the data is obtained using a pure 
bremsstrahlung model with a plasma temperature of 14.1^3 q keV consistent, within the large 
errors, with the 9.6±0.3keV value measured by BeppoSAX. 
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In addition, we analyzed a 10 ks XMM- Newton observation together with the BAT data. 
Utilizing the surface brightness profile obtained by lO Tiara et all (120041 ) (beta model with 
/?=0.55 and core radius R c =2.6' ) we determine that our region of 10' radius includes ~90 % of 
the cluster emission. We employ such cross-normalization factor when fitting XMM- Newton 
and BAT data. The BAT-XMM-Newton spectra, shown in the right panel of Fig. [T3J, are 
well fitted by an APEC model with a plasma temperature of 9.27lg;27keV and abundance 
of 0.25(±0.04) solar. The 99% upper limit on the 2-10 keV non-thermal flux of 2.70 xl0~ 12 
erg cm~ 2 s _1 is in disagreement with t he non-thermal flux of (4 .0±0.1)xl0~ 11 erg cm" 2 s _1 



detected in the same band by RXTE (IGruber fc Rephaelill2002l ) 



Harris & Mileyl (119781 ) disco vered a diffuse radio halo associated with the A 2319 cluster. 
An intensive study was done by lFeretti et all (Il997l ). from which we take S 61Q mh z = Uy 



and a = 0.92 to estimate the lower limit on the magnetic field reported in Table [3j 



3. Clusters Properties 

3.1. Constraints on non-thermal excess emission 

In order to constrain the non-thermal hard X-ray emission, we have produced 3 a upper 
limits on the 50-100 keV non-thermal flux for each source presented in the previous section. 
We excluded the Perseus and the Coma clusters. Indeed, Perseus is the only cluster where 
the detected "hard-tail" is certainly produced by the brightest AGN while Coma requires a 
dedicated analysis. We chose the 50-100 keV energy band because above 50keV the thermal 
emission of the clusters is negligible. 

The 3er upper limit has been computed integrating the source flux in the 50-100 keV 
range and subtracting the thermal flux arising from the best thermal fit. We added to this 
value three times the 1 a uncertainty. The upper limits are reported in Table [3J These upper 
limits were derived using BAT data alone. It is important to note that, indeed, thanks 
to the very good sensitivity of BAT, all these upper limits are very stringent. Indeed, the 
non-thermal flux for all these sources is constrained to be below ~1 mCrab. 

In the derivation above, we do not make any assumption on the mechanism generating 
the non-thermal flux. Hi 
pal emission proce ss (e.g. 



Sarazin 


1999: 


Nevalainen et al. 


2004 



2007 



Eckert et al.ll2008l ). If this is true, then the IC emission ca n be modeled as a power-law 
with photon index ~2 in the 1-200 keV energy range (see e.g. iReimer et al.ll2004l ). We thus 
computed the 99 % CL upper limits to the IC flux in the 50-100 keV band by adding a 
power-law model to the best fits reported in Tab. [2j These limits are reported in Tab. HI It 
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is worth noting that, since we are using XMM-Newton/XRT/Chandra and BAT data, these 
upper limits are a factor 5-10 lower than those derived using BAT data alone (see Tab. [3]). 



3.2. Stacking analysis 

A few clusters show positive, marginal, residuals above 50 keV; this is the case for 
Abell 3266, Abell 3571 and Abell 2142. Such features are not statistically significant to 
warrant an additional component (e.g. non-thermal power-law). However, it might be that 
the non-thermal component be just below the BAT sensitivity for such clusters. In this 
case the stacking technique offers the capability to explore the average properties of a given 
population beyond the current instrumental limit. Thus, we produced the stacked spectrum 
of all clusters except Perseus and Coma (for the reasons explained above). The average 
spectrum is produced by the weighted average of all the spectra. The weight is chosen to be 
the inverse of the variance of a given bin and it is exactly the same procedure used to extract 
the spectra of each individual source. The same stacking t echnique has been applied with 



success to the study of Seyfert galaxies detected by BAT (lAjello et al.l l2008d ). The total 
spectrum has an exposure time of ~56 Ms and it is shown in Fig. [TH A fit with a simple 
bremsstrahlung model yields a good chi-square (x 2 /dof = 7.2/10). The best fit temperature 
is 10.8^0 8 keV which is in very good agreement with the mean temperature of 10.4 keV as 
determined by averaging the values obtained fitting a simple bremsstrahlung model to each 
cluster's spectrum (using BAT data alone). This is a good confirmation that the chosen 
stacking technique reproduces well the average properties of our cluster sample. 

From the best thermal fit, we derive a 99 % CL upper limit (50-100 keV) for the 
non-thermal component of 1.9xl0~ 12 erg cm" 2 s" 1 (0.3mCrab). At the average redshift 
of the sample (z=0.05 8), this translates into a limiting luminosity of 1.4xl0 43 erg s _1 . 
Nevalainen et al. J2004h reported the detection of an average non-thermal component de- 



tected in the stacked spectrum (20-80 keV) of BeppoSAX clusters. Their non-thermal lu- 
minosity is comprised □ in the (0.5-5.0) xlO 43 erg s" 1 range. In the 20-80 keV band our 
99 % CL limit on the non-thermal luminosity is 2.2xl0 43 erg s _1 . Thus, the findings of 



Nevalainen et al.l (I2004J ) are consistent with our analysis. 



On the other hand, all clusters, except perhaps Perseus and Abell 3571, are undergoing 
a merging phase. These last two clusters are those which show the lowest ICM temper- 
ature in our sample. The L x — T relation (shown in Fig. [T5l) reinforces the picture that 



3 The measurement reported by 
in this paper. 



Nevalainen et al. 



(120041 ) had to be converted to the Hubble constant used 
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most of the BAT clusters are mergers. Indeed, the best fit to the data with a power-law 
of the form L = AqTq, where T 6 = T/ 6keV (fixin g a at 2.8£ yields a norma l izatio n 



A 6 =(2.82±0.8)xl0 43 h^erg s" 1 . While. Markevitchl (119981 ) and lArnaud fc Evrardl f|1999h . 



found for A 6 a value of (12.53±1.08)xl0 43 erg s" 1 and (12.13±0.06) xlO 43 erg s" 1 
respectively. Indeed, merging cluste rs are known to s egregate at lower luminosities (or higher 
temperatures) in the L x — T plane (jOta et al.ll2006l ). 



There is a growing evidence which points to wards a rather non-uniform distribution of 



see 



Markevitch et al. 


2003; 


O'Hara et al. 


2004 



Eckert et al.ll20070. Both hy drodynamical simulations (e.g. 



Takizawa 



19991 ) and observations 



Markevitch et al.ll2003l . for Abell 0754) have shown that shocks due to cluster mergers 
can heat the ICM up to ~15 keV. Figure [TH] shows that for the merging clusters the mean 
temperature measured by BAT is slightly larger (given the large uncertainties) than the 
mean ICM temperature measured below lOkeV. A similar trend, althoug h using different 



wave bands, has been recently reported for a sample of 192 galaxy clusters (jCavagnolo et al. 



20081 ). Moreover, for the merging clusters, the BAT centroid is shifted to positions where 
Chandra and XMM- Newton have detected regions of hot gas. Based on these evidences, we 
believe that the conjecture that these clusters show regions of "hot" gas is a more viable 
claim than the one which foresees the presence of a strong IC component. 

This claim is also supported by the fact that the high-energy residuals (e.g. residuals 
above 10 keV of the spectral fits using a single thermal model) are in general better described 
by an additional thermal component than a power-law model. To prove this, we selected 
those clusters which show, in the analysis presented in § 12.21 the largest residuals above 10 keV 
from the thermal model used. These clusters, which are Abell 2029, Triangulum Australis 
and Abell 2319, also show a large deviation between the ICM temperature measured below 
and above lOkeV (see Fig. [JB]). We made a fit to each of these clusters with: 1) a single 
thermal model, 2) the sum of a thermal and a power-law model, and 3) the sum of two 
thermal models. The residuals to each of these fits are shown in Figures [T71 [TBI an d HH] 
while the spectral parameters are summarized in Table O We note that in all three cases the 
additional thermal model explains the residuals better than an additional power-law model. 
We also remark that for most of the BAT clusters (in this case for Triangulum Australis and 
Abell 2319) the single thermal model is already a good description of the data (x 2 /dof=~1.0) 
and given the statistics no other additional model is required. This means that currently 
the high-energy residuals (with respect a single thermal fit) are not significant. Longer BAT 
exposures will clarify the existence and nature of these emissions. 



4 Given the small ran ge in luminosity spanned by our sample we fixed a at the value determined by 
Arnaud fc Evrardl j 19991 ). 



-22 - 



3.3. Cluster Magnetic Field Assessment 



The diffuse synchrotron radio emission (radio halos, relics and mini-halos) proves the 
existence of magnetic fields in the ICM. The intensity of the synchrotron emission depends 
both on the strength of the magnetic field and on the electron density. If the non-thermal 
X-ray emission results from IC scattering off the same radio electrons by the CMB, then the 
degen eracy in magnetic field and relativistic electron density can be broken (e.g. iRephaeli 
200ll ) . Therefore the non-detection of a non-thermal component can be used to place a lower 
limit on the magneti c fields B in clusters (the rati o of I C to radio flux is inversely proportional 
to B a+1 ). Following lHarris fc Romanishinl (119741 ) and ISarazinl (119881 ). we estimate the lower 
limit on B (the volume averaged component along the line of sight): 



fx 



2.47 x 10 



19 n MB Kp) 



r 

J v, 



B a(p) 



4960 T CMB 
B 



where a is the spectral index, p = 2 a + 1, f x the X-ray flux integrated over the band 



between v min = 50 KeV and v max = 100 KeV (f x = k c f"' 



*di/ x , in erg cm 2 s *) 



k s v 01 the flux density at the radio frequency u, 



m erg cm s 



Hz" 



T. 



CMB 



2.7 'K 



the temperature of the CMB and a(p) and b(p) as in ISarazinl (119881 ) (eq. 5.6 and 5.8). Since 
our clusters are nearby, in the formula above we neglect redshift corrections. 

Although the limit on the X-ray flux is very stringent, the measurement of the diffuse 
radio emission is complicated by the presence of individual radio galaxies in the cluster. In 
most cases the radio observations were not sensitive enough over a wide range of spatial 
scales to subtract the contribution of the single sources. Moreover, the spectral index varies 
with the distance from the center of the cluster. These factors make the derivation of the 
magnetic field intensity uncertain. Therefore the values listed in table 3 have to be taken 
as order of magnitude estimates. Such estimates point to magnetic fields that are typically 
a fraction of a fiG. These low values indicate that these systems are far from equipartition. 
This is possible if one considers that the magnetic fields and the relativistic particles may 
have a different spatial extension and history. 

The magnetic field can also be evaluated by measuring the Faraday rotation (RM) of 



Kim et al. 


1991; 


Clarke et al. 


2001) 



likely because the interpretation of Faraday rotation measurements and the derivation of 
the mean magnetic field strength rests on assumptions of the magnetic field topology (see 



Goldshmidt &: Rephaelilll993l ; IColafrancesco et al 



20051 . for an extensive discussion) . We can 



produce a more robust, upper limit on the IC flux considering that the IC emission spectrum 
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can be approximated as a power law in the 1-200 keV energy band (see e.g. iReimer et al. 



2004 . for more details). Using both 2-10 keV and BAT data we were able to produce the 



limits reported in Table H] which are in some cases a factor 5-10 lower than our previous esti- 
mated values (Tab. [3]) based on BAT data alone. This in turn translates in larger intensities 
of the magnetic field which in a few cases reach the ~ 0.5 //G value. 



4. Cluster Population 
4.1. Cluster log N - log S 



Thanks to the serendipitous character of the BAT survey, it is possible to derive, for 
the first time, the source counts distribution (also known as log N - log S) of clusters above 
15 keV. This can be obtained as: 



8=0 1 



(2) 



where N$ is the total number of clusters with fluxes greater than S and f2j is the geometrical 
area surveyed to that limiting flux. The cumulative distribution is reported in Fig.[20j Source 
counts distributions are generally fitted by a power law of the form N(> S) = AS~ a . Given 
the small number of objects, we do not attempt a maximum likelihood fit to derive the slope 
a, but we note that our flux distribution is consistent with an Euclidean function N oc S~ 3 ' 2 
as shown in Fig. [201 We derive the normalization A as that one which reproduces the number 



of observed objects above the flux o f ~ 1 x 1 11 e rg cm 2 s 1 . Using the 90% confidence 



limits for small numbers derived by iGehrelsl (119861 ). we find that a good representation of 
our data is obtained by iV(> S) = (4. 19^4 x 10~ 4 deg _2 )S' ] ~ 1 1 ' 5 where Sn is the flux in unit 
of 10 -11 erg cm -2 s -1 . This function is also shown in Fig. | 



Interestingly, we note that the integrated flux of all clusters above 10 -11 erg cm -2 s _1 is 
9.7xl0~ n erg cm -2 s _1 sr _1 . This is only ~0.1% o f the Cosmic X-ray Background (CXB) 
flux as measured by BAT in the 15-5 5 keV band (lAiello et al.l l2008bl ). but 5-10% of the 



total flux resolved by BAT into AGNs ( Ajello et al. 2008cl ). Thus, clusters of galaxies 



are a 



sizable population among the extragalactic objects (mostly AGN) detected by BAT. 



We can compare the BAT log N - log S, with those derived in the 0.5-2 keV band. 
In doing so, we extrapolate the BAT spectra to the 0.5-2 keV band using the temperatures 
measured below 10 keV. The cluster surface density above 10 -12 erg cm -2 s _1 in the 0.5- 
2.0 keV band is 4.3jjj x 1 0~ 2 deg" 2 whic h is in rather good agreement with the finding of 



Vikhlinin et all (11998 ) and iBurenin et al.l (120071 ). 
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The BAT source counts distribution can be used to estimate the foreseen number of 
galaxy clusters above a given flux limit. In doin g so, we adopt for a the -1.4 value which has 



been established by deeper X-ray surveys (e.g. IJones et al 



19981 : IBohringer et alJl200lL and 



references therein). Indeed, using the -3/2 value would certainly overestimate the cluster 
density at lower fluxes. As an example, an instrument surveying the whole sky to 10~ 13 erg 
cm -2 s _1 would detect approximately ~10000 galaxy clusters in the 15-55 keV band. The 
BAT sample itself will comprise up to 30 objects, if BAT will able to reach the 0.5mCrab 
flux limit on the whole sky. 



4.2. X— ray Luminosity Function 

Since all our clusters have a measured redshift, we can derive their luminosity function. 
Its construction relies on the knowledge of the survey volume V max as a function of X-ray 
luminosity. The survey volume is the volume of the cone defined by the survey area and the 
luminosity distance at which a cluster with a given luminosity could just be observed at the 
flux limit. The limiting luminosity distance Dl u m , and thus also V max , can be determined 
solving iteratively the following equation: 

LUm AnF lim k(T,z) {) 

where L x is the source luminosity and k(T, z) is the ^-correction which accounts for the 
redshift ing of the source spectrum. 

Once the V max is computed for each object, the cumulative luminosity function can be 
derived as: 

- 1 

N(>L X ) = J2 77 — TTT Mo Mpc- 3 ] . (4) 
The cumulative luminosity fu nction of the BAT clus ters, obtained with the method 



reported above, is shown in Fig. [211 IBohringer et al.l (120021 ). analyzing a flux-limited sample 



of ROSAT galaxy clusters (REFLEX), derived that a good parametrization of the differential 
luminosity function is a Schechter function of the form: 

— = n exp(--) (-)'"—. (5) 
In order to compare the REFLEX l uminosity function with the BAT one, we adopt for no, 



and a, the values determined by IBohringer et al.l (120021 ). Moreover, since the REFLEX 



luminosity function is derived in the 0.1-2.4 keV band, we need to convert the luminosities to 
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the BAT, 15-55 keV, band. We do this using the mean clusters temperature (kT = 8.1keV) 
determined in the 2-10 keV band (see right panel of Fig. [T6|) . The reason for adopting this 
temperature instead of the BAT-derived temperature is double. First, given the S/N ratio, 
temperatures determined in the 2-10 keV band have a better accuracy than temperatures 
determined in the BAT band. Most importantly however, using the 2-10 keV temperature 
allows a more accurate extrapolation of the source luminosity from the ROSAT (0.1-2.4 keV) 
to the BAT (15-55 keV) band. The extrapolated, cumulative, REFLEX luminosity function 
is also reported Fig. [2TJ It is apparent that, notwithstanding the extrapolation, the agree- 
ment of the BAT data and the REFLEX luminosity function is excellent. This agreement 
is not however surprising because most of the BAT clusters constitute the bright end of 
the REFLEX luminosity function. The value of converted to the 15-55 keV band is 
L,=7.3xl0 43 h^erg s" 1 while n = S.13±H x 10~ 7 and a = 1.63 

Integrating the luminosity function multiplied by the luminosity yields the total X- 
ray emissivity W of galaxy clusters. Above the survey limit of 2 x 10 43 erg s" 1 , we fmc|f| 
1^=2.83 xlO 37 erg s _1 Mpc~ 3 ( 15-55 keV). This can be compared to the total emissivity of 
AGN which w a s der ived for the local Universe and a similar energy band (17-60 keV) by 



Sazonov et al.l (120071 ). After correcting for the small difference between the energy bands, 
the AGN local emissivity above 2 x 10 43 erg s -1 is M/^ G jv=14.1 xlO 37 erg s -1 Mpc~ 3 . It is 
thus clear that galaxy clusters contribute substantially (~20 % level with respect to AGN) 
to the local X-ray output. 



Discussion 



5.1. Non-thermal hard X-ray emission 



Direct evidence of the presence of relativistic electrons in the IC M arises from the 
existence of large radio halos (jDennison!ll980l ; iFeretti fc Giovanninill2007l ) . The same electron 
population responsible for the synchrotron emission can in principle scatter CMB photons 
by IC and produce hard X-ray radiation. The intensity of this radiation relative to the 
synchrotron emission ultimately depends on the value of the magnetic field. 

A firm detection of non-thermal components in the spectra of galaxy clusters has re- 
mained elusive in the past as well as in this study. Indeed, Perseus is the only galaxy cluster 
in the BAT sample where a non-thermal high-energy component is revealed at high signifi- 
cance. Most likely this component is due to the emission of the central AGN NGC 1275. The 



We do not provide an error estimate since the luminosity function was not fitted to the data. 
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rest of the clusters detected by BAT do not show a significant non-thermal emission. Using 
BAT data alone, we are able to constrain the non-thermal component below the mCrab level 
in the 50-100 keV energy band. The BeppoSAX detection above 50 keV of an average non- 
thermal co mponent in the stacked spectrum of several clusters is consistent with the BAT 
upper limit (INevalainen et al.ll2004f). As discussed i n § 12.21 some of the individual detections 
of non-thermal components (e.g. lEckert et al.ll2008l ) are consistent (albeit some marginally) 
with the upper limits derived using BAT data alone. Thus, we cannot exclude that such 
non-thermal components exist and that they are currently below or at the limit of the BAT 
sensitivity. If we assume that the principal emission mechanism is IC scattering of GeV elec- 
trons off CMB photons, then the cluster magnetic field is constrained to be >0.1 /iG. These 
low magnetic intensities would show that the magnetic field is far from equipartition (i.e. 
the ene rgy in the magnetic field is different with respect to the electrons energy). As pointed 
out by iPetrosian fc Bykovl (120081 ). this can happen if the sources generating the magnetic 
field and accelerating the electrons are not identical. 

However, IC emission by relativistic electrons can be modeled as a p ower-law in the 



1-200 keV energy regime (e.g. iNevalainen et al.l 120041 ; iReimer et al.l 120041 ) . Thus, using 



XMM-Newton/XRT/Chandra and BAT data we are able to constrain, more robustly, the 
IC emission mechanism. With this approach, we confirm the detection and the flux of the 
hard component in the spectrum of Abell 0754, but we are also able to prove (thanks to 
the resolution of XMM-Newton) that a single point-like object, 2MASS 09091372-0943047, 
located less than 2' from the BAT centroid accounts for the whole non-thermal emission. 
For the rest of the clusters, we are able to produce upper limits which are a factor 5- 
10 lower than previously estimated. These limits in turn translate into a slightly larger 



Kim et al. 


1991; 


Clarke et al. 


2001) 



then the chances of d etecting IC emission from clusters with the currently flying instruments 



become really small (IPfrommer fc Enfilinl 120041 ) . Indeed, the values of the predicted IC flux 



account for only <10 % of the claimed non-thermal X-ray emission abov e lOkeV when takin g 



both pri mary and seconda ry-generated electrons into account (see e.g. iMiniati et al.ll200ll ). 



Recently, IPfrommer! (120081 ). using high- resolution simulations of a sample of representative 
galaxy clusters, showed that the predicted IC flux for the Coma and Perseus clusters would 
be a factor 50 lower than the claimed detections. 

Our combined analysis, thus, put tight constrains on the IC mechanism. However, IC 
emission is the process that most likely explains the claimed non-thermal emission, but not 
the only one. Hard X-ray flux from galaxy clusters can be interpreted as bremsstrahlung from 



acceleration in the turbulent ICM (e.g. 



Enfilin et al. 


1999; 


Petrosian 


2001) 
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the radio and the non-thermal X-ray flux are no longer strictly related and equipartition may 
apply. However, the non-thermal bremsstrahlung model requires a continuos input of energy 
in the ICM which as a consequence will cause its tempe rature to increase. Thus, the non- 
thermal bremsstrahlung phase is likely to be short lived (IPetrosian fc Bykovll2008l ). 



5.2. Structure Formation 



All the galaxy clusters detected by BAT, except perhaps Abell 3571, are merging sys- 
tems. Some, as Abell 0754, Abell 2142 and, Abell 3266, are experiencing violent merging due 
to encounters of subclusters w ith comparable ma s ses. In the commo n scenarios of hierar- 
chical structure formation (e.g. iMiniati et al.ll2000i ; iRyu fc Kang 120031). large system s evolve 
as the result of merging of smaller structures. As reviewed in iDolag et al.l (120081 ). cluster 
mergers generate internal shocks ( Mach number less than 4) which provide most of the ICM 
gas heating (e.g. iQuilis et al.lll998l ). and also likely convert a non-negligible fraction (<10 %) 
of their power into CRs. The shocks primarily heat the ions because the kinetic en ergy of an 
ion e ntering the shock region is larger than that o f an electron by their m ass ratio (ITakizawa 
19991 ). Cosmological simulations have shown (e.g. iPfrommer et al.l 120071 ). that in the case of 



ongoing merger activity, the relative CR pressure (to the thermal ICM pressure) is greatly 
enhanced, up to 15-20 %, due to strong merger shock waves. This pressure is likely larger 
in the outskirts of the cluster because of the lower sound speed and t he larger density of th e 
ICM in the central region which makes CR acceleration less efficient ( IPfrommer et al.ll2007l ). 



Hot spots as well as cold fronts hav e been found in many merging clusters thanks 
to the superior resolution of Chandra (e.g. iMarkevitch et al.ll2000l ; iMarkevitch fc Vikhlinin 
200ll ; IMarkevitch et al.ll2003l ). Hydrodynamical simulations have highlighted that ~1 Gyr 
after the encounter of two clusters with compar able masses po s t-shock regions with hig h 



temperatures (T ?«10-20keV ) are formed (e.g. iTakizawal Il999l ; iRitchie &: Thomasl 120021 ). 
In the BAT sample there is a clear correlation of gas temperature and merging activity. 
Indeed, Abell 3571 and Perseus, which are in a late merging stage, display the lowest plasma 
temperatures among the clusters in our sample. INTEGRAL recently unveiled the presence 
of a hotter region (T= 12 ± 2 keV) located south-west of the centre of the Coma cluster 
( lEckert et al.l 120071 ) . These findings highlight the important role of merging shocks in the 
heating of ICM. 
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5.3. Clusters Statistics 



The serendipitous character of the BAT survey allowed us to determine, for the first 
time above lOkeV, the log N- log S and luminosity function distributions of galaxy clusters. 
Both are in very good agreement with previous studies. The log N - log S highlights that 
the clusters BAT detects produce a negligible fraction (~0.1%) of the X-ray background 
emission, but they represent a sizable population (5-10 %) with respect to the local AGN. 
The BAT log N- log S shows that future instruments with a sensitivity 10 or 100 times better 
than BAT (above 15keV) will detect clusters at a density of ~0.01deg -2 and ~0.24deg -2 
respectively. 

The BAT luminosity distribution allowed us to determine that the volume emissivity of 
galaxy clust ers is W(> 2 x 10 43 er g s~ 1 )=2.38xl0 37 erg s _1 Mpc~ 3 . Above the same limiting 



luminosity, ISazonov et al. J2007h derived that the volume emissivity of the local AGN is 
W^4G7v=14.1xl0 37 erg s _1 Mpc -3 . Thus, above 2 x 10 43 erg s _1 , the cluster volume emissivity 
is 20 % of that one of AGN. Integrating the luminosity functions to lower luminosity (e.g. 
10 41 erg s -1 ) changes this fraction to ~10 %. This change is due to the fact that at low 
luminosity the AGN luminosity function is steeper than the cluster luminosity function (e.g. 
Sazonov et all 120071 : iBohringer et alJbooj ). 



5.4. Future Prospects 



The study of non-thermal processes in clusters of galaxies requires a multi-wavelength 
approach. The ongoing Swift/BAT survey will likely comprise up to 30 clusters if an all-sky 
sensitivity of 0.5mCrab will be reached and it will improve the signal-to-noise ratio for the 
spectra of the clusters presented here. Ultimately, major progresses are expected with the 
launch of Simbol-X@, XEUS0, NUSTAR0, and NeXlfl. Indeed, their sensitivities and 
spectro-imaging capabilities up to high energies (80 keV and beyond) will provide new and 
better constraints on the hard X-ray emission. 

The future generation of radio arrays combined with high-energy observations will allow 
to shed some light on the energetics of relativistic particles, the nature and frequency of 



6 http://www. asdc.asi.it/simbol-x/ 

7 http: / / www.rssd.esa.int /index. php?project=xeus 

8 http:/ /www. nustar.caltech.edu/ 

9 http:/ /www. astro. isas.ac.jp/future/NeXT/ 
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acceleration processes, and the strength and structure of magnetic fields. As we already 
discussed, this astrophysical information has strong cosmological implications. The Long 
Wavelength ArraJ^I (LWA), the Low Frequency Array!"] (LOFAR), and ultimately the 
Square Kilometre ArrajEl (SKA), will operat e over a critical radio frequency range to detect 
relativistic plasma in large-scale structure and clusters in a sensitive way. The advance in 
sensitivity and resolution will increase the statistics of known radio halos and radio relics at 
different redshif ts. The correlation of sensit ive X-ray and radio detections will be particularly 
important (e.g. lEnfilin fc Rottgering!l2002l ). At the same time, thanks to the high angular 
and spectral resolution, the Faraday rotation studies will significantly improve yielding a 
better determination of the cluster magnetic field. 

Much attention is directed towards the Gamma-ray Large Area Space Telescope @(GLAST) 
which, with an unprecedented sensitivity, spatial resolution and dynamic range at GeV en- 
ergies, will shed light on the origin of the extragalactic 7-ray background. Galaxy clusters 
and shocks from st ructure format ions are natural candidates for explaining part of this dif- 
fuse emission (e.g. iDermerl 120071 . and references therein). All the BAT clusters are good 
candidate for GLAST since they are nearby and are mergers. Indeed, in merg ing systems. 



part of the internal shocks energy is very likely co nverted into CRs acceleration (IDolag et al. 



20081 ). As pointed out by iPfrommer et al.l (120081 ). above 100 MeV the cluster emission will 
likely be dominated by pion decay 7-rays even though a contribution from non-thermal 
bremsstrahlung and IC emission of secondary electrons is expected. This will provide a 
unique information about the hadron component of CRs which is not included in estimates 
of CR pressure based only on the observations discussed above concerning electrons and 
magnetic field. Since cosmic ray protons loss time is long, the 7r°-bump detection would 
proof that hadrons ha ve been confined in the ICM for as long as the Hubble time (e.g. 



Berezinsky et al.lll997l ). Stringent constraints on the cosmic rays content in the ICM is 



fundamental for the future space missions which will use galaxy clusters to constrain and 
understand the nature of Dark Energy (e.g. eROSITA0). 



10 http: / /lwa. unm.edu 

n http: / / www.lofar.org/ 

12 http:/ /www. skatelescope.org 

13 http:/ /www-glast. stanford.edu 

1 http://www.mpe. mpg.de/projects. html^crosita 
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6. Conclusions 

BAT is the first instrument to detect above 15 keV an all-sky sample of galaxy clusters^. 
The BAT energy range (15-200 keV) is the best one to investigate the presence of non-thermal 
emission whose detection remained so far controversial. The results of our investigation can 
be summarized as follows: 



• Perseus is the only cluster among the 10 BAT objects which displays an high-energy 
non-thermal component which extends up to 200 keV. It is very likely that the cen- 
tral AGN NGC 1275 is responsible for such emission. Thi s claim is supported b y 



several evidences: 1) the variability se en with BeppoSAX ([Nevalainen et al.l 12004 ) 



2) the XMM-Newton spectral analysis (IChurazov et al.l 120031 ) . and 3) our combined 



BAT-XRT-XMM-Newton analysis which shows that the nucleus has a typical AGN 
spectrum. 

• The BAT spectra of the remaining 9 galaxy clusters is well fitted by a simple thermal 
model that constrains non-thermal flux to be below 1 mCrab in the 50-100 keV band. 

• Assuming that IC scattering is the main mechanism at work for producing non-thermal 
high-energy flux, it is possible to estimate the magnetic field using Radio data and the 
upper limits derived above. We obtain that in all the BAT clusters the (average) 
magnetic field is > 0.1 fiG. These (rather uncertain) values are in disagreement (if the 
magnatic field intensities are close to the lower limits) with the, also uncertain, Faraday 
rotation measurements which show that the magnetic field is in the ~ fiG range. Our 
low magnetic field values would imply that the magnetic field is far from equipartition. 

• The stacked spectrum of the BAT clusters (except Perseus and Coma) confirms once 
again the absence of any non-thermal high-energy component. The ~56Ms stacked 
spectrum constrains any non-thermal flux to be below 0.3mCrab (or 1.9xl0 _12 erg 
cm~ 2 s _1 ) in the 50-100 keV band. 

• Using Swift/XRT, XMM-Newton and Chandra, in addition to BAT data, we were able 
to produce X-ray cluster spectra which extend more than 3 decades in energy (0.5- 
50keV). In all cases, but Perseus and Abell 0754, the broad-band X-ray spectrum is 
well approximated by a single-temperature thermal model. These spectra allowed us 
to put constrains on the IC emission mechanism which are a factor >5 lower than 



15 We are aware of an independent work (jOkajima et al.1 120081 ) based on an alternative analysis of BAT 
survey data which reaches conclusions consistent with this analysis. 
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those derived using BAT data alone. This would in turn imply a larger intensity of the 
magnetic field. For both Perseus and Abell 0754 an additional power-law component is 
statistically required, but several evidences confirm that two X-ray point sources (NGC 
1275 and 2MASS 09091372-0943047) account for the total non-thermal emission. 

• The cluster centroid shift in different wavebands, the morphology and the complex 
temperature maps (available in literature), show that 8 out of 10 clusters are in the 
middle of a major merging phase. Shocks, which are revealed by XMM-Newton and 
Chandra images, are actively heating the ICM as the BAT high temperatures testify. 
The BAT observations and limits on the non-thermal emissions can help to calibrate 
the large scale structure formation simulations focusing in particular on the treatment 
of non-thermal particle emission and cooling. 

• We have produced the first cluster source count (also known as log N - log S) dis- 
tribution above 15keV. This shows that, at the limiting fluxes sampled by BAT, the 
surface density of clusters is ~5 % of that one of AGNs. Moreover, we find that the 
contribution of clusters to the Cosmic X-ray background is of the ~0.1 % order in the 
15-55 keV band. The BAT log N - log S can be used to predict the cluster surface 
density for future hard X-ray instruments. 

• The X-ray luminosity function of the BAT clusters, the first derived above 15keV, is 
in excellent agreement with the ROSAT luminosity function derived in the 0. 1-2.4 keV 
band. 



MA acknowledges funding from the DFG Leibniz-Prize (HA 1850/28-1). PR is sup- 
ported by the Pappalardo Postdoctoral Fellowship in Physics at MIT. NC was partially 
supported from a NASA grant NNX07AV03G. MA and PR wish to acknowledge Balu for 
his incomparable constant enthusiasm. We thank T. Okajima for providing a copy of his 
manuscript before submission and for interesting discussions. The anonymous referee is 
aknowledged for his helpful comments which improved the manuscript. This research has 
made use of the NASA/IPAC extragalactic Database (NED) which is operated by the Jet 
Propulsion Laboratory, of data obtained from the High Energy Astrophysics Science Archive 
Research Center (HEASARC) provided by NASA's Goddard Space Flight Center, of the 
SIMBAD Astronomical Database which is operated by the Centre de Donnees astronomiques 
de Strasbourg, and of the ROSAT All Sky Survey maintained by the Max Planck Institut 
fur Extraterrestrische Physik. 



Facilities: Swift (BAT/XRT), XMM-Newton, Chandra 



-32 - 



REFERENCES 

Ajello, M., Greiner, J., Kanbach, G., Rau, A., Strong, A. W., & Kennea, J. A. 2008a, ApJ, 
678, 102 

Ajello, M., Greiner, J., Sato, G., & Kanbach, G. 2008b, ApJ, submitted 

Ajello, M., Rau, A., Greiner, J., Kanbach, G., Salvato, M., Strong, A. W., Barthelmy, S. D., 
Gehrels, N., Markwardt, C. B., & Tueller, J. 2008c, ApJ, 673, 96 

Ajello, M., Rebusco, P., Cappelluti, N., Greiner, J., & Reimer, O. 2008d, (in preparation) 

Anders, E. & Grevesse, N. 1989, Geochim. Cosmochim. Acta, 53, 197 

Ando, S. & Nagai, D. 2007, ArXiv e-prints, 705 

Arnaud, M., Aghanim, N., Gastaud, R., Neumann, D. M., Lumb, D., Briel, U., Altieri, B., 
Ghizzardi, S., Mittaz, J., Sasseen, T. P., k Vestrand, W. T. 2001, A&A, 365, L67 

Arnaud, M. & Evrard, A. E. 1999, MNRAS, 305, 631 

Avni, Y. 1976, ApJ, 210, 642 

Bacchi, M., Feretti, L., Giovannini, G., & Govoni, F. 2003, A&A, 400, 465 

Barthelmy, S. D., Barbier, L. M., Cummings, J. R., Fenimore, E. E., Gehrels, N., Hullinger, 
D., Krimm, H. A., Markwardt, C. B., Palmer, D. M., Parsons, A., Sato, G., Suzuki, 
M., Takahashi, T., Tashiro, M., & Tueller, J. 2005, Space Science Reviews, 120, 143 

Berezinsky, V. S., Blasi, P., & Ptuskin, V. S. 1997, ApJ, 487, 529 

Bird, A. J., Barlow, E. J., Bassani, L., Bazzano, A., Belanger, G., Bodaghee, A., Capitanio, 
F., Dean, A. J., Fiocchi, M., Hill, A. B., Lebrun, F., Malizia, A., Mas-Hesse, J. M., 
Molina, M., Moran, L., Renaud, M., Sguera, V., Shaw, S. E., Stephen, J. B., Terrier, 
R., Ubertini, P., Walter, R., Willis, D. R., & Winkler, C. 2006, ApJ, 636, 765 

Blanton, E. L., Sarazin, C. L., & McNamara, B. R. 2003, ApJ, 585, 227 

Blasi, P., Gabici, S., & Brunetti, G. 2007, ArXiv Astrophysics e-prints 

Bohringer, H., Collins, C A., Guzzo, L., Schuecker, P., Voges, W., Neumann, D. M., 
Schindler, S., Chincarini, G., De Grandi, S., Cruddace, R. G., Edge, A. C, Reiprich, 
T. H., & Shaver, P. 2002, ApJ, 566, 93 



-33- 

Bohringer, H., Schuecker, P., Guzzo, L., Collins, C. A., Voges, W., Schindler, S., Neumann, 
D. M., Cruddace, R. G., De Grandi, S., Chincarini, G., Edge, A. C, MacGillivray, 
H. T., & Shaver, P. 2001, A&A, 369, 826 

Bonamente, M., Lieu, R., & Mittaz, J. P. D. 2001, ApJ, 561, L63 

Bowyer, S., Berghdfer, T. W., & Korpela, E. J. 1999, ApJ, 526, 592 

Brown, D. L. & Burns, J. O. 1991, AJ, 102, 1917 

Brunetti, G., Blasi, P., Cassano, R., & Gabici, S. 2004, MNRAS, 350, 1174 

Burenin, R. A., Vikhlinin, A., Hornstrup, A., Ebeling, H., Quintana, H., & Mescheryakov, 
A. 2007, ApJS, 172, 561 

Bykov, A. M., Bloemen, H., & Uvarov, Y. A. 2000, A&A, 362, 886 

Catura, R. G, Fisher, P. G, Johnson, H. M., & Meyerott, A. J. 1972, ApJ, 177, LI 

Cavagnolo, K. W., Donahue, M., Voit, JVL, & Sun, M. 2008, ApJ(accepted) 

Chen, Y., Reiprich, T. H., Bdhringer, H., Ikebe, Y., & Zhang, Y.-Y. 2007, A&A, 466, 805 

Churazov, E., Forman, W., Jones, G, & Bdhringer, H. 2003, ApJ, 590, 225 

Churazov, E., Forman, W., Vikhlinin, A., Tremaine, S., Gerhard, O., & Jones, C. 2007, 
ArXiv e-prints, 711 

Clarke, T. E., Blanton, E. L., & Sarazin, C. L. 2004, ApJ, 616, 178 

Clarke, T. E., Kronberg, P. P., & B6hringer, H. 2001, ApJ, 547, Llll 

Colafrancesco, S., Marchegiani, P., & Perola, G. C. 2005, A&A, 443, 1 

Condon, J. J., Cotton, W. D., Greisen, E. W., Yin, Q. F., Perley, R. A., Taylor, G. B., & 
Broderick, J. J. 1998, AJ, 115, 1693 

Condon, J. J., Griffith, M. R., & Wright, A. E. 1993, AJ, 106, 1095 

Dar, A. & Shaviv, N. J. 1995, Physical Review Letters, 75, 3052 

De Grandi, S. & Molendi, S. 1999, ApJ, 527, L25 

-. 2002, ApJ, 567, 163 



-34- 



de Plaa, J., Werner, N., Bykov, A. M., Kaastra, J. S., Mendez, M., Vink, J., Bleeker, 
J. A. M., Bonamente, M., & Peterson, J. R. 2006, A&A, 452, 397 

Dennison, B. 1980, ApJ, 239, L93 

Dermer, C. D. 2007, ApJ, 659, 958 

Dolag, K., Bykov, A. M., & Diaferio, A. 2008, Space Science Reviews, 29 

Durret, F., Slezak, E., Lieu, R., Dos Santos, S., & Bonamente, M. 2002, A&A, 390, 397 

Ebisawa, K., Bourban, G., Bodaghee, A., Mowlavi, N., & Courvoisier, T. J.-L. 2003, A&A, 
411, L59 

Eckert, D., Produit, N., Neronov, A., & Courvoisier, T. J. . 2007, ArXiv e-prints, 710 

Eckert, D., Produit, N., Paltani, S., Neronov, A., & Courvoisier, T. J.-L. 2008, A&A, 479, 
27 

Edge, A. C, Stewart, G. C, & Fabian, A. C. 1992, MNRAS, 258, 177 

Ensslin, T. A., Biermann, P. L., Kronberg, P. P., & Wu, X.-P. 1997, ApJ, 477, 560 

Enfilin, T. A., Lieu, R., & Biermann, P. L. 1999, A&A, 344, 409 

Enfilin, T. A. & Rdttgering, H. 2002, A&A, 396, 83 

Ettori, S., Fabian, A. C, & White, D. A. 1998, MNRAS, 300, 837 

Fabian, A. C, Pringle, J. E., & Rees, M. J. 1976, Nature, 263, 301 

Felten, J. E., Gould, R. J., Stein, W. A., & Woolf, N. J. 1966, ApJ, 146, 955 

Feretti, L. & Giovannini, G. 1998, in Untangling Coma Berenices: A New Vision of an Old 
Cluster, ed. A. Mazure, F. Casoli, F. Durret, & D. Gerbal, 123-+ 

Feretti, L. & Giovannini, G. 2007, preprint (astro-ph/0703494) 

Feretti, L., Giovannini, G., & Bohringer, H. 1997, New Astronomy, 2, 501 

Finoguenov, A., Briel, U. G., Henry, J. P., Gavazzi, G., Iglesias-Paramo, J., & Boselli, A. 
2004, A&A, 419, 47 

Finoguenov, A., Henriksen, M. J., Miniati, F., Briel, U. G., & Jones, C. 2006, ApJ, 643, 790 
Forman, W., Kellogg, E., Gursky, H., Tananbaum, H., & Giacconi, R. 1972, ApJ, 178, 309 



-35 - 



Fritz, G., Davidsen, A., Meekins, J. F., & Friedman, H. 1971, ApJ, 164, L81 

Fujita, Y., Hayashida, K., Nagai, M., Inoue, S., Matsumoto, H., Okabe, N., Reiprich, T. H., 
Sarazin, C. L., & Takizawa, M. 2008, ArXiv e-prints, 806 

Fujita, Y., Kohri, K., Yamazaki, R., k Kino, M. 2007, ApJ, 663, L61 

Fujita, Y., Takizawa, M., & Sarazin, C. L. 2003, ApJ, 584, 190 

Fusco-Femiano, R., dal Fiume, D., Feretti, L., Giovannini, G., Grandi, P., Matt, G., Molendi, 
S., & Santangelo, A. 1999, ApJ, 513, L21 

Fusco-Femiano, R., Landi, R., & Orlandini, M. 2007, ApJ, 654, L9 

Fusco-Femiano, R., Orlandini, M., Brunetti, G., Feretti, L., Giovannini, G., Grandi, P., & 
Setti, G. 2004, ApJ, 602, L73 

Fusco-Femiano, R., Orlandini, M., De Grandi, S., Molendi, S., Feretti, L., Giovannini, G., 
Bacchi, M., & Govoni, F. 2003, A&A, 398, 441 

Gehrels, N. 1986, ApJ, 303, 336 

Gehrels, N., Chincarini, G., Giommi, P., Mason, K. O., Nousek, J. A., Wells, A. A., White, 
N. E., Barthelmy, S. D., Burrows, D. N., Cominsky, L. R., Hurley, K. C, Marshall, 
F. E., Meszaros, P., Roming, P. W. A., Angelini, L., Barbier, L. M., Belloni, T., 
Campana, S., Caraveo, P. A., Chester, M. M., Citterio, O., Cline, T. L., Cropper, 
M. S., Cummings, J. R., Dean, A. J., Feigelson, E. D., Fenimore, E. E., Frail, D. A., 
Fruchter, A. S., Garmire, G. P., Gendreau, K., Ghisellini, G., Greiner, J., Hill, J. E., 
Hunsberger, S. D., Krimm, H. A., Kulkarni, S. R., Kumar, P., Lebrun, F., Lloyd- 
Ronning, N. M., Markwardt, C. B., Mattson, B. J., Mushotzky, R. F., Norris, J. P., 
Osborne, J., Paczynski, B., Palmer, D. M., Park, H.-S., Parsons, A. M., Paul, J., 
Rees, M. J., Reynolds, C. S., Rhoads, J. E., Sasseen, T. P., Schaefer, B. E., Short, 
A. T., Smale, A. P., Smith, I. A., Stella, L., Tagliaferri, G., Takahashi, T., Tashiro, 
M., Townsley, L. K., Tueller, J., Turner, M. J. L., Vietri, M., Voges, W., Ward, M. J., 
Willingale, R., Zerbi, F. M., & Zhang, W. W. 2004, ApJ, 611, 1005 

Giacconi, R., Murray, S., Gursky, H., Kellogg, E., Schreier, E., & Tananbaum, H. 1972, 
ApJ, 178, 281 

Giovannini, G. & Feretti, L. 2000, New Astronomy, 5, 335 

Giovannini, G., Feretti, L., Venturi, T., Kim, K.-T., & Kronberg, P. P. 1993, ApJ, 406, 399 



-36- 



Gisler, G. R. & Miley, G. K. 1979, A&A, 76, 109 

Gitti, M., Brunetti, G., & Setti, G. 2002, A&A, 386, 456 

Goldshmidt, O. & Rephaeli, Y. 1993, ApJ, 411, 518 

Gruber, D. & Rephaeli, Y. 2002, ApJ, 565, 877 

Harris, D. E., Bahcall, N. A., & Strom, R. G. 1977, A&A, 60, 27 

Harris, D. E. & Miley, G. K. 1978, A&AS, 34, 117 

Harris, D. E. & Romanishin, W. 1974, ApJ, 188, 209 

Henriksen, M. J. & Markevitch, M. L. 1996, ApJ, 466, L79 

Henry, J. P. & Briel, U. G. 1995, ApJ, 443, L9 

-. 1996, ApJ, 472, 137 

Henry, J. P., Finoguenov, A., & Briel, U. G. 2004, ApJ, 615, 181 

Ikebe, Y., Reiprich, T. H., Bdhringer, H., Tanaka, Y., & Kitayama, T. 2002, A&A, 383, 773 
Inoue, S., Aharonian, F. A., & Sugiyama, N. 2005, ApJ, 628, L9 

Johnston, M. D., Bradt, H. V., Doxsey, R. E., Marshall, F. E., Schwartz, D. A., & Margon, 
B. 1981, ApJ, 245, 799 

Jones, B. B. & Finlay, E. A. 1974, Australian Journal of Physics, 27, 687 

Jones, L. R., Scharf, G, Ebeling, H., Perlman, E., Wegner, G., Malkan, M., & Horner, D. 
1998, ApJ, 495, 100 

Kaastra, J. S., Bykov, A. M., Schindler, S., Bleeker, J. A. M., Borgani, S., Diaferio, A., 
Dolag, K., Durret, F., Nevalainen, J., Ohashi, T., Paerels, F. B. S., Petrosian, V., 
Rephaeli, Y., Richter, P., Schaye, J., & Werner, N. 2008, Space Science Reviews, 36 

Katz, J. I. 1976, ApJ, 207, 25 

Kempner, J. C. & Sarazin, C. L. 2001, ApJ, 548, 639 

Kennea, J. A., Markwardt, C. B., Tueller, J., Barthelmy, S., Burrows, D., Campana, S., 
Falcone, A., Gehrels, N., Grupe, D., Godet O., Krimm, H., Mangano, V., Mineo, T., 
Pagani, G, Page, K., Mushotzky, R. F., & Skinner, G. K. 2005, The Astronomer's 
Telegram, 677, 1 



-37- 



Kim, K.-T., Kronberg, P. P., k Tribble, P. C. 1991, ApJ, 379, 80 

Lieu, R., Mittaz, J. P. D., Bowyer, S., Lockman, F. J., Hwang, C.-Y., k Schmitt, J. H. M. M. 
1996, ApJ, 458, L5 

Lumb, D. H., Warwick, R. S., Page, M., & De Luca, A. 2002, A&A, 389, 93 

Lutovinov, A. A., Vikhlinin, A., Churazov, E. M., Revnivtsev, M. G., & Sunyaev, R. A. 
2008, ArXiv e-prints, 802 

Markevitch, M. 1998, ApJ, 504, 27 

Markevitch, M., Forman, W. R., Sarazin, C. L., & Vikhlinin, A. 1998, ApJ, 503, 77 

Markevitch, M., Mazzotta, P., Vikhlinin, A., Burke, D., Butt, Y., David, L., Donnelly, H., 
Forman, W. R., Harris, D., Kim, D.-W., Virani, S., & Vrtilek, J. 2003, ApJ, 586, L19 

Markevitch, M., Ponman, T. J., Nulsen, P. E. J., Bautz, M. W., Burke, D. J., David, L. P., 
Davis, D., Donnelly, R. H., Forman, W. R., Jones, C, Kaastra, J., Kellogg, E., Kim, 
D.-W., Kolodziejczak, J., Mazzotta, P., Pagliaro, A., Patel, S., Van Speybroeck, L., 
Vikhlinin, A., Vrtilek, J., Wise, M., & Zhao, P. 2000, ApJ, 541, 542 

Markevitch, M. & Vikhlinin, A. 2001, ApJ, 563, 95 

Markevitch, M. L., Sarazin, C. L., & Irwin, J. A. 1996, ApJ, 472, L17 

Markwardt, C. B., Tueller, J., Skinner, G. K., Gehrels, N., Barthelmy, S. D., & Mushotzky, 
R. F. 2005, ApJ, 633, L77 

Mills, B. Y., Slee, O. B., & Hill, E. R. 1960, Australian Journal of Physics, 13, 676 

Miniati, F., Jones, T. W., Kang, H., k Ryu, D. 2001, ApJ, 562, 233 

Miniati, F., Ryu, D., Kang, H., Jones, T. W., Cen, R., k Ostriker, J. P. 2000, ApJ, 542, 608 
Miralda-Escude, J. k Babul, A. 1995, ApJ, 449, 18 
Molendi, S. k De Grandi, S. 1999, A&A, 351, L41 

Molendi, S., de Grandi, S., Fusco-Femiano, R., Colafrancesco, S., Fiore, F., Nesci, R., k 
Tamburelli, F. 1999, ApJ, 525, L73 

Molendi, S. k Gastaldello, F. 2008, arXiv:0807.2653 

Nagai, D., Vikhlinin, A., k Kravtsov, A. V. 2007, ApJ, 655, 98 



-38- 



Nevalainen, J., Kaastra, J., Parmar, A. N., Markevitch, M., Oosterbroek, T., Colafrancesco, 
S., & Mazzotta, P. 2001, A&A, 369, 459 

Nevalainen, J., Oosterbroek, T., Bonamente, M., & Colafrancesco, S. 2004, ApJ, 608, 166 

O'Hara, T. B., Mohr, J. J., & Guerrero, M. A. 2004, ApJ, 604, 604 

Okajima, T., Tueller, J., & Gehrels, N. 2008, (in preparation) 

Ota, N., Kitayama, T., Masai, K., & Mitsuda, K. 2006, ApJ, 640, 673 

Peres, C. B., Fabian, A. C, Edge, A. C, Allen, S. W., Johnstone, R. M., & White, D. A. 
1998, MNRAS, 298, 416 

Perkins, J. S., Badran, H. M., Blaylock, G., Bradbury, S. M., Cogan, P., Chow, Y. C. K., Cui, 
W., Daniel, M. K., Falcone, A. D., Fegan, S. J., Finley, J. P., Fortin, P., Fortson, L. F., 
Gillanders, G. H., Gutierrez, K. J., Grube, J., Hall, J., Hanna, D., Holder, J., Horan, 
D., Hughes, S. B., Humensky, T. B., Kenny, G. E., Kertzman, M., Kieda, D. B., 
Kildea, J., Kosack, K., Krawczynski, H., Krennrich, F., Lang, M. J., LeBohec, S., 
Maier, G., Moriarty, P., Ong, R. A., Pohl, M., Ragan, K., Rebillot, P. F., Sembroski, 
G. H., Steele, D., Swordy, S. P., Valcarcel, L., Vassiliev, V. V., Wakely, S. P., Weekes, 
T. C, & Williams, D. A. 2006, ApJ, 644, 148 

Petrosian, V. 2001, ApJ, 557, 560 

Petrosian, V. & Bykov, A. M. 2008, Space Science Reviews, 25 
Pfrommer, C. 2008, MNRAS, 385, 1242 
Pfrommer, C. & EnBlin, T. A. 2004, A&A, 413, 17 
Pfrommer, C, EnBlin, T. A., & Springel, V. 2008, MNRAS, 277 

Pfrommer, C, EnBlin, T. A., Springel, V., Jubelgas, M., & Dolag, K. 2007, MNRAS, 378, 
385 

Primini, F. A., Basinska, E., Howe, S. K., Lang, F., Levine, A. M., Lewin, W. H. G., 
Rothschild, R., Baity, W. A., Gruber, D. E., Knight, F. K., Matteson, J. L., Lea, 
S. M., & Reichert, G. A. 1981, ApJ, 243, L13 

Quilis, V., Ibanez, J. M. A., & Saez, D. 1998, ApJ, 502, 518 

Read, A. M. & Ponman, T. J. 2003, A&A, 409, 395 



-39- 

Reimer, A., Reimer, O., Schlickeiser, R., & Iyudin, A. 2004, A&A, 424, 773 

Reimer, O., Pohl, M., Sreekumar, P., & Mattox, J. R. 2003, ApJ, 588, 155 

Renaud, M., Belanger, G., Paul, J., Lebrun, F., & Terrier, R. 2006a, A&A, 453, L5 

Renaud, M., Gros, A., Lebrun, F., Terrier, R., Goldwurm, A., Reynolds, S., & Kalemci, E. 
2006b, A&A, 456, 389 

Rephaeli, Y. 1979, ApJ, 227, 364 

Rephaeli, Y. 2001, in American Institute of Physics Conference Series, Vol. 558, American 
Institute of Physics Conference Series, ed. F. A. Aharonian & H. J. Volk, 427 

Rephaeli, Y. & Gruber, D. 2002, ApJ, 579, 587 

Rephaeli, Y., Nevalainen, J., Ohashi, T., & Bykov, A. M. 2008, Space Science Reviews, 134, 
71 

Revnivtsev, M., Sazonov, S., Jahoda, K., & Gilfanov, M. 2004, A&A, 418, 927 
Ritchie, B. W. & Thomas, P. A. 2002, MNRAS, 329, 675 
Robertson, J. G. & Roach, G. J. 1990, MNRAS, 247, 387 
Rossetti, M. & Molendi, S. 2004, A&A, 414, L41 
-. 2007, preprint (astro-ph/0702417) 

Ryu, D. & Kang, H. 2003, Journal of Korean Astronomical Society, 36, 105 
Sanders, J. S. & Fabian, A. C. 2007, MNRAS, 381, 1381 

Sanders, J. S., Fabian, A. C, Allen, S. W., & Schmidt, R. W. 2004, MNRAS, 349, 952 

Sanders, J. S., Fabian, A. C, & Dunn, R. J. H. 2005, MNRAS, 360, 133 

Sanderson, A. J. R., Ponman, T. J., & O'Sullivan, E. 2006, MNRAS, 372, 1496 

Sarazin, C. L. 1988, X-ray emission from clusters of galaxies (Cambridge Astrophysics Series, 
Cambridge: Cambridge University Press, 1988) 

-. 1999, ApJ, 520, 529 

Sarazin, C. L. & Kempner, J. C. 2000, ApJ, 533, 73 



-40- 



Sarazin, C. L., Wise, M. W., & Markevitch, M. L. 1998a, ApJ, 498, 606 
-. 1998b, ApJ, 498, 606 

Sauvageot, J. L., Belsole, E., & Pratt, G. W. 2005, A&A, 444, 673 

Sazonov, S., Revnivtsev, M., Krivonos, R., Churazov, E., & Sunyaev, R. 2007, A&A, 462, 
57 

Slee, O. B. 1977, Australian Journal of Physics Astrophysical Supplement, 43, 1 

Slee, O. B. & Higgins, C. S. 1975, Australian Journal of Physics Astrophysical Supplement, 
36, 1 

Smith, R. K., Brickhouse, N. S., Liedahl, D. A., & Raymond, J. C. 2001, ApJ, 556, L91 

Stephen, J. B., Bassani, L., Malizia, A., Bazzano, A., Ubertini, P., Bird, A. J., Dean, A. J., 
Lebrun, F., & Walter, R. 2006, A&A, 445, 869 

Takizawa, M. 1999, ApJ, 520, 514 

Taylor, G. B., Barton, E. J., & Ge, J. 1994, AJ, 107, 1942 

Thierbach, M., Klein, U., & Wielebinski, R. 2003, A&A, 397, 53 

Timokhin, A. N., Aharonian, F. A., & Neronov, A. Y. 2004, A&A, 417, 391 

Tueller, J., Barthelmy, S., Burrows, D., Falcone, A., Gehrels, N., Grupe, D., Kennea, J., 
Markwardt, C. B., Mushotzky, R. F., & Skinner, G. K. 2005a, The Astronomer's 
Telegram, 668, 1 

-. 2005b, The Astronomer's Telegram, 669, 1 

Valinia, A., Henriksen, M. J., Loewenstein, M., Roettiger, K., Mushotzky, R. F., & Madejski, 
G. 1999, ApJ, 515, 42 

Venturi, T., Bardelli, S., Zagaria, M., Prandoni, I., & Morganti, R. 2002, A&A, 385, 39 

Vikhlinin, A., McNamara, B. R., Forman, W., Jones, C, Quintana, H., & Hornstrup, A. 
1998, ApJ, 498, L21 

Voges, W., Aschenbach, B., Boiler, T., Brauninger, H., Briel, U., Burkert, W., Dennerl, 
K., Englhauser, J., Gruber, R., Haberl, F., Hartner, G., Hasinger, G., Kurster, M., 
Pfeffermann, E., Pietsch, W., Predehl, P., Rosso, C, Schmitt, J. H. M. M., Triimper, 
J., & Zimmermann, H. U. 1999, A&A, 349, 389 



-41 - 

Watanabe, M., Yamashita, K., Furuzawa, A., Kunieda, H., & Tawara, Y. 2001a, PASJ, 53, 
605 

-. 2001b, PASJ, 53, 605 

Werner, N., Kaastra, J. S., Takei, Y., Lieu, R., Vink, J., & Tamura, T. 2007, A&A, 468, 849 
White, III, R. E., Day, C. S. R., Hatsukade, I., & Hughes, J. P. 1994, ApJ, 433, 583 
Willson, M. A. G. 1970, MNRAS, 151, 1 



This preprint was prepared with the AAS IATgX macros v5.2. 



-42- 




Fig. 1. — Assessment of systematic errors for two representative energy channels: 18-22 keV 
(left) and 57.6-75.4 keV (right). The histograms show the distribution of S/N for 160 random 
positions (noise) in the sky away from known or detected sources. The dashed line is a fit 
to the data using a Gaussian profile. The 1 a widths of the Gaussian profiles are compatible 
with 1.0. 
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Fig. 2. — 3(T average spectral sensitivity as a function of energy based on the analysis of 160 
randomly extracted spectra. The dashed line is the Crab Nebula spectrum divided by 1000. 
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Fig. 3. — Le/£ Panel: ROSAT 0. 1-2.4 keV surface brightness of Perseus with BAT signif- 
icance contours superimposed. The contours range from 2.5 a to 28 a. Right Panel: Joint 
XRT-BAT spectrum of Perseus. The best fit (thick solid line) is the sum of two thermal 
models (dashed and long-dashed line) and of a power-law component (thin solid line). 
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Fig. 4.- Pane/: ROSAT 0.1-2.4 keV surface brightness of Abell 3266 with BAT 

significance contours superimposed. The contours range from 2.5cr to 5.5 a. Right Panel: 
Joint fit to XMM-Newton-BAT data for Abell 3266 with a thermal model. The best model 
is shown as solid line. 
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Fig. 5.- Left Panel: ROSAT 0.1-2.4 keV surface brightness of Abell 0754 with BAT 
significance contours superimposed. The contours range from 2.5cr to 8.0cr. Right Panel: 
Joint fit to XMM-Newton-BAT data. The best fit model (thick solid line) is the sum of a 
thermal model (dashed line) and of a power law (thin solid line). 




Fig. 6. — Left Panel: ROSAT 0. 1-2.4 keV surface brightness of Coma with BAT significance 
contours superimposed. The contours range from 2.5 a to 20 a. Right Panel: Joint fit to 
XMM-Newton-BAT data. The best fit model (solid line) is the sum of two thermal models 
(dashed and dotted lines). 
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Fig. 7. — Residuals to the fit to Coma data using: a single thermal model (top), sum of a 
thermal model and a power law (middle), and the sum of two thermal models (bottom). 
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Fig. 8.- Left Panel: ROSAT 0.1-2.4 keV surface brightness of Abell 3571 with BAT 
significance contours superimposed. The contours range from 2.5 cr to 5.0 a. Right Panel: 
Joint fit to XMM-Newton-BAT data with a thermal model. The best fit model is shown as 
a solid line. 
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Fig. 9.- Pane/: ROSAT 0.1-2.4 keV surface brightness of Abell 2029 with BAT 

significance contours superimposed. The contours range from 2.5 cr to 5.0 cr. Right Panel: 
Joint XRT-BAT spectrum of Abell 2029. The best fit (thick solid line) is the sum of two 
thermal models (thin solid and dashed line). 



-47- 




Fig. 10.- Left Panel: ROSAT 0. 1-2.4 keV surface brightness of Abell 2142 with BAT 
significance contours superimposed. The contours range from 2.5 o to 7.0 a. Right Panel: 
Joint fit to XMM-Newton-BAT data for Abell 2142 with single thermal model. The best fit 
model is shown as a solid line. 




Fig. 11. — Left Panel: XMM-Newton 1.0-7.0 keV surface brightness of the Triangulum Aus- 
tralis cluster with BAT significance contours superimposed. The contours range from 2.5 ex 
to 7.0 a. Right Panel: Joint fit to XMM-Newton-BAT data for the Triangulum Australis 
cluster with a thermal model. The best fit model is shown as a solid line. 
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Fig. 12. — Left Panel: ROSAT 0.1-2.4 keV surface brightness of the Ophiucus cluster with 
BAT significance contours superimposed. The contours range from 2.5 o to 22 a. Right 
Panel: Joint fit to Chandra-BAT data for the Ophiucus cluster with a thermal model. The 
best fit is shown as solid line. 




Fig. 13.- Left Panel: ROSAT 0. 1-2.4 keV surface brightness of Abell 2319 with BAT 
significance contours superimposed. The contours range from 2.5 o to 22 o. Right Panel 
Joint fit to XMM-Newton-BAT data for Abell 2319 Australis cluster. The best fit model 
thermal model is shown as a solid line. 
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Fig. 14. — Stacked spectrum of the clusters in our sample and the best fit (dashed line) 
with a bremsstrahlung model. 
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Fig. 15. — Luminosity- Temperature relation for the BAT cluste rs. The black line is the 
best, power-law, fit to the data while the gray line is the best fit of lArnaud fc Evrardl (jl999j) 
converted to the BAT energy band. 
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Fig. 16. — Comparison of best fit gas temperatures. The x-axis reports the temperatures 
derived using BAT data (above 15keV) while the y-axis shows the temperatures derived 
using 2-10 keV data (XMM-Newton, Chandra or XRT). The dashed line shows the Tbat — 
^2-iofceV function. The largest deviations are for the merging clusters Abell 2029 and Abell 
2319. 
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Fig. 17. — Residuals to the fit to Abell 2029 data using: a single thermal model (top), sum 
of a thermal model and a power law (middle), and the sum of two thermal models (bottom). 
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Fig. 18. — Residuals to the fit to Triangulum Australis data using: a single thermal model 
(top), sum of a thermal model and a power law (middle), and the sum of two thermal models 
(bottom). 
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Fig. 19. — Residuals to the fit to Abell 2319 data using: a single thermal model (top), sum 
of a thermal model and a power law (middle), and the sum of two thermal models (bottom). 
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Fig. 20. — Cumulative flux number relation for the BAT clusters (15-55 keV). The dashed 
line is an overlaid power law N(>S) = A S~ L5 . 




Fig. 21. — Cumulative luminosity function of the BAT clusters (15- 55 keV). The solid line 
is the X-ray luminosity function determined for the REFLEX survey (jBohringer et al.ll2002l ) 
converted to the BAT energy band. 
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Table 1. Clusters detected in the 15-55 keV band 



NAME 


R.A. 


DEC 


S/N 


ID 


z 


EXPOSURE 


OFFSET 




(J2000) 


(J2000) 








(Ms) 


(arcmin) 


J0319.8+4130 


49.9573 


41.5110 


28.00 


Perseus 


0.0175 


2.89 


0.5 


J0431.3-6126 


67.8297 


-61.4388 


5.61 


Abell 3266 


0.0590 


3.81 


2.1 


J0908.9-0938 


137.2391 


-9.6346 


8.28 


Abell 0754 


0.0530 


2.96 


1.8 


J1259.4+2757 


194.8531 


27.9523 


19.95 


Coma Cluster 


0.0230 


4.32 


5.1 


J1347.7-3253 


206.9500 


-32.9000 


5.05 


Abell 3571 


0.0397 


1.78 


4.5 


J1511.0+0544 


227.7500 


5.7485 


5.33 


Abell 2029 


0.0770 


2.71 


0.8 


J1558.5+2714 


239.6256 


27.2417 


7.11 


Abell 2142 


0.0890 


3.62 


3.3 


J1638.8-6424 


249.7136 


-64.4000 


6.90 


Triangulum A. 


0.0510 


1.77 


4.9 


J1712.3-2319 


258.0914 


-23.3242 


21.63 


Ophiucus 


0.028 


1.30 


1.7 


J1920.9+4357 


290.2405 


43.9646 


11.72 


Abell 2319 


0.056 


3.87 


2.2 
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Table 2. Spectral parameters from combined XMM-Newton/XRT/Chandra and BAT fits 

(errors are 90% C.L.) 



NAME Flux a L x a kT T model X 2 /dof 

(10" 11 cgs) (10 43 crg s" 1 ) (keV) 



Perseus 


3.90^° 


2.7±?;J 


3.00+^?/6.40 + °- 6 7 ? 


1.71q;7 apec+apec+pow 


152.8/144 


Abell 3266 


0.73jft2 


6 q +0 ' 9 

u - y -0.9 




apec 


666.8/841 


Abell 0754 


i 1 1 +0.04 
± - ii -0.04 


0-^-0.3 


q. 3+0-4 


apec+pow 


1217.0/1072 


Coma b 


o OO+0.23 
Z.OO_ 22 


3 +0 - 2 

O.U_q 4 


8.40±«;^7i a4 45^;- 


apec+apec 


846.5/856 


Abell 3571 


fiS +a09 

u -°°-0.06 


9 7+0.3 
z - ' -0.4 




apec 


723.9/1367 


Abell 2029 


1 m+ - 16 


16.8^ 


4.1±^/9.6±I;g 


apec+apec 


394.2/363 


Abell 2142 


o qn+ 010 

u - yu -0.10 


21-5±|;| 


8-40182 


apec 


361.9/398 


Triangulum A. 


i - ou -0.10 


o O+0.6 
O-O-0.2 


q qn+ - 30 
y - ou -0.30 


apec 


925.8/1074 


Ophiucus 


c 7+0.5 
°-'-0.5 


O QQ+0.28 
y - oo -0.14 


q qq+o.24 
y.»o_ 24 


apec 


323.1/351 


Abell 2319 


i,oo -0.14 


1 3 +0 - 9 


9 23+ - 27 

J.ZO„ 27 


apec 


1151.3/1274 



a Flux and Luminosities are computed in the 15-55 keV band. 

b The spectral values reported for Coma are only representative for the source extraction region (i.e. 
10' around the BAT centroid; see § 12.4.41 for more details). 
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Table 3. 3<x Upper limits on the non-thermal component and Clusters' properties 



NAME 


CC a ? Merger? F 50 _ioo^ b 


B 


^radio ^radio 


a 


Ref 




(10~ 12 erg cm -2 s^ 1 ) 




(Jy) (MHz) 







Perseus 


y 


y 














Abell 3266 


n 


y 


<5.30 


> 0.17 


1.070 


2700 


0.95 


1 


Abell 0754 


n 


y 


<6.50 


> 0.10 


0.086 


1365 


1.5 


2 


Coma 


n 


y 














Abell 3571 


y d 


n e 


<11.5 


> 0.03 


0.0084 


1380 


1.5 9 


3 


Abell 2029 


j d 


y 


<4.83 


> 0.25 


0.528 


1380 


1.5 9 


4 


Abell 2142 


y f 


y 


<5.35 


>0.06 


0.0183 


1400 


1.5 9 


5 


Triangulum A. 


y f 


y / 


<4.65 


> 0.17 


< 0.033 


4850 


1.5 9 


6 


Ophiucus 


n 


n 


<5.89 


> 0.11 


6.4 


160 


2.0 


7 


Abell 2319 




y 


<3.41 


> 0.10 


1.0 


610 


0.92 


8 



a CC=Cool Core 

b BAT data alone were used to estimate the upper limits 
c References for the radio flux 
d Moderate CC 

e The morphology and temperature map indicate that it is a relaxed cluster, but the radio structure 
points at late stages of merging 

f Under discussion 

g Arbitrary spectral index 
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Table 4. Non-thermal emission from combined XMM-Newton/XRT/Chandra and BAT 

data. 



NAME 


F50-100 keV a 


B b 




(10~ 12 erg cm 2 s" 1 ) 


( M G) 


Perseus 






Abell 3266 


<0.57 


> 0.55 


Abell 0754 






Coma 






Abell 3571 


1 4 +o.4 


~ 0.08 


Abell 2029 


<1.27 


> 0.42 


Abell 2142 


<1.50 


> 0.10 


Triangulum Australis 


<0.65 


> 0.39 


Ophiucus 


<2.80 


> 0.15 


Abell 2319 


<0.67 


> 0.15 



a The flux has been estimated using a power-law 
spectrum with a photon index of 2.0 in the 1- 
200 keV energy band. Upper limits are 99 % CL 
while errors are 90 % CL. 

b In order to compute the intensity of the mag- 
netic field we used the same Radio data reported 
in Tab. [3] 
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Table 5: Comparison of different spectral fits to the clusters which show a large deviation 
between the ICM temperature as measured below and above 10 keV. As a reference for the 
reader also the parameters of Coma are reported. kTi and kT 2 are the temperatures of the 
two thermal models (in keV) while norm, and T are the normalization at 1 keV (in ph cm -2 
s _1 keV -1 ) and the photon index of the power-law model. Frozen parameters do not have 
an error estimate. 



Cluster 


Thermal 


Thermal + power law 


Thermal + Thermal 


Abell 2029 








kT, 


6 7^+0-52 
°- ' °-0.31 


6 78+ - 46 

°- ' °-0.33 


4t -- L — 1.5 


r 




2.0 




norm. 




1.55+H? x 10~ 3 

— 1.15 




kT 2 






9.6 +2 /° 

^ w — 2.0 


X 2 /dof 


407.3/364 


402.1/363 


394.2/363 


Triangulum A. 








kT, 


q Qn+0-30 


y --^ u -0.28 


ii 1+0.34 
ii - i -0.27 


r 




2.0 




norm. 




< 1.40 x 10~ 4 




kT 2 






l.OO_ 27 


X 2 /dof 


925.8/1074 


925.8/1073 


916.5/1072 


Abell 2319 








kT x 

r 

norm. 


O 90+0.27 
y - z,o -0.27 


Q QQ+0.35 
a -°°-0.52 

i 7+0.2 

l - ' -0.3 

7.81^ x 10~ 4 


1 1 2 +as 

il - z -1.0 


kT 2 






1 q+i-64 

L -^-0A0 


X 2 /dof 


1151.34/1274 


1139.9/1272 


1127.8/1272 


Coma 








kT x 

r 

norm 


6 sn +0 - 09 

°- ou -0.05 


7 iq+0.16 
' - ly -0.06 
o i 1 +0.10 
z - 11 -0.13 

3.56±°;|t x 10~ 3 


8.4o±8:lt 


kT 2 






1.45lg;S 


X 2 /dof 


1168.9/858 


905.5/856 


846.5/856 



